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I.  INTRODUCTION 


Various  physical  and  engineering  problems  require  the  evaluation  of 
integrals,  over  the  real  line,  of  complex  functions  that  have  singularities, 
such  as  poles  and  branch  points,  along  the  path  of  integration.  These 
integrals,  which  appear  frequently  when  solving  differential  or  integral 
equations  by  integral  transform  methods,  ordinarily  arise  when  the  inverse 
transforms  need  to  be  evaluated. 

As  is  well-known,  such  integrals  may  have  no  meaning  at  all  in  the  sense 
of  Riemann's  theory  of  integration.  The  frustrated  user  of  the  integral 
transform  method  may  at  this  point  become  discouraged,  if  not  despairing, 
and  forsake  this  approach  for  more  conventional,  though  less  direct,  methods. 

This  state  of  affairs  is  common  to  other  areas  of  human  activity  and  may 
well  stem  from  the  implicit  faith  one  has  in  the  power  of  a particular 
discipline  to  produce  automatically  the  correct  results  when  applied  to  the 
problems  of  a different  discipline.  Mathematics  in  itself,  when 
legitimately  applied,  will  yield  not  only  the  correct  result,  but  also  the 
entire  class  of  all  possible  results  to  the  problem  under  consideration. 

This  clearly  implies  that  we  cannot  live  by  mathematics  alone.  Our 
salvation  lies  then  in  choosing  that  solution,  from  the  set  of  all 
mathematically  possible  ones,  which  satisfies  the  basic  principles  of  the 
discipline  which  gave  rise  to  the  original  problem.  To  illustrate,  it  is 
well-known  that  waves  cannot  propagate  inside  a quiescent  medium  prior  to 
the  activation  of  a source  of  disturbances.  Nevertheless,  one  of  the 
mathematical  solutions  #to  the  wave  equation  appropriate  to  this  problem  is  a 
wave  propagating  before  the  disturbances  have  been  initiated.  This  wave, 
known  as  the  advanced  solution,  is  usually  discarded,  since  it  does  not 
conform  with  our  perception  of  physical  causality,  i.e.,  that  effect  should 
not  precede  cause.  Examples  abound  in  the  physical  and  engineering 
literature  where  certain  possible  solutions  are  discarded  because  these 
violate  basic  requirements  of  the  various  theories  to  which  the  problems 
belong. 

In  the  present  report,  the  author  proposes  to  list  the  possible 
mathematical  interpretations  of  the  integrals  mentioned  above,  restricted  at 
present  only  to  poles  along  the  path  of  integration,  none  of  which  occur  at 
an  endpoint  of  the  interval  of  integration. 

We  start  this  report  in  Section  II  with  a discussion  of  integrals  having 
simple  poles  at  a finite  distance  from  the  origin,  the  interval  of 
integration  being  the  whole  real  axis.  We  shall  call  these  full-range 
integrals.  In  three  subsections  we  proceed  to  give  three  interpretations, 
that  increase  in  generality,  as  follows:  The  Principal  Value  of  Cauchy,  the 
Plemelj  formulas  obtained  by  path  deformation,  and  the  6-function 
interpretation  which  includes  the  first  two  as  special  cases.  It  is  of 
interest  to  note  that  the  third  interpretation  has  not  been  sufficiently 
emphasized  in  the  literature,  although  it  is  based  on  well-known  results  of 
generalized  function  theory  (the  solution  of  equations  with  vanishing 
coefficients,  or  the  zero  divisor  problem). 
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In  Section  III  we  consider  full-range  integrals  possessing  multiple 
poles  on  the  real  axis.  The  results  of  Section  II  are  extended,  with  the 
appropriate  interpretation,  also  to  this  case.  The  need  for  a more  general 
definition  of  the  Principal  Value  is  shown  by  exhibiting  the  failure  of  the 
conventional  one  to  exist  for  a simple  example. 

Similar  results  are  presented  in  Sections  IV  and  V for  half-range 
integrals  and  finite  intervals.  In  the  treatment  are  included  both  simple 
and  multiple  poles.  In  these  sections  several  methods  are  proposed  for  the 
explicit  evaluation  of  the  Principal  Value,  while  these  and  other  methods 
are  described  in  Appendix  I for  the  full  range  integrals.  Some  of  these 
methods  may  be  novel.  The  fact  of  the  matter  is  that  the  literature  on  the 
subject  is  rather  silent  about  effecting  the  evaluation  of  the  Principal 
Value  by  methods  which  avoid  the  e-limit  process.  As  can  be  seen  from  the 
text,  the  Principal  Value,  when  it  exists,  is  made  the  cornerstone  of  all 
the  possible  interpretations  of  a singular  integral  and  thus  it  is  essential 
to  have  an  explicit  expression  for  it,  whether  in  closed  form  or  as  an 
ordinary,  convergent  integral. 

In  Appendix  II  the  theory  is  applied  to  specific  examples  for  full-range 
integrals  with  simple  poles.  We  mention  in  particular  the  results  obtained 
when  there  is  an  infinite  number  of  poles.  The  method  used  is  of  interest 
in  itself,  being  essentially  an  extension,  to  certain  classes  of  functions, 
of  the  partial  fraction  decomposition  formulas  for  polynomials  (Ref.  10). 
Since  infinite  sums  occur  in  this  case,  questions  of  convergence  are 
raised.  It  is  not  surprising,  perhaps,  that  these  can  be  bypassed  by  an 
appeal  to  the  theory  of  generalized  functions.  This  theory  also  provides 
the  interpretation  when  the  full-range  integrals  do  not  converge,  even  with 
the  poles  removed.  Thus,  although  the  Principal  Value  does  not  exist  in  the 
conventional  sense,  it  does  so  as  a generalized  function.  This  aspect  of 
the  Principal  Value  does  not  seem  to  have  been  highlighted  in  previous  works. 

In  Appendix  III  one  example  of  a full-range  integral  with  a multiple 
pole  is  presented.  Appendices  IV  and  V treat  examples  for  half-range  and 
finite  intervals,  respectively,  with  simple  poles.  Examples  for  multiple 
poles  have  not  been  given  for  these  cases,  but  their  treatment  presents  no 
particular  difficulty. 

In  the  concluding  remarks,  a short  outline  of  the  main  theme  of  the  work 
is  given.  Several  related,  important  topics,  not  investigated  here,  which 
require  a similar  systematic  treatment,  are  mentioned. 


II.  INTEGRALS  WITH  SIMPLE  POLES 
In  this  section  we  shall  consider  integrals  of  the  type 

w,-  , n _ r<*»  f(x.X) 


X)  - C 
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In  eq.  (1),  f(x,X)  is  a complex  function  of  the  reel  variables  x and  X. 
The  latter  can  be  viewed  as  a varying  paraneter,  independent  of  x.  In 
addition, £ is  real  and  f(x,X)  is  defined  for  x*5,  with  f(5, X)  ■ 0.  Hence 
the  integrand  in  I(5,X)  has  a real  pole  of  first  order  at  x*5. 

I(5,X),  as  it  stands,  represents  three  limit  processes,  each  of  which 
may  yield  a divergent  result.  From  Analysis  we  know  that  eq.  (1)  is  a 
shorthand  notation  for 

1(5, X)  - lim  Jb  dx  (2) 

a-*— oo  a *”  ^ 
b-*-*-ao 

where  a -*-•  and  b independently.  The  integral  over  the  finite  interval 
a,b  represents  the  third  limit  process  (one  of  summation  in  this  case) 
alluded  to  above.  Let  us  denote  this  integral  by  l]j.  Due  to  the 
presence  of  the  pole  at  x*5,  this  integral  does  not  exist  in  the 
conventional  sense. 


In  the  following  we  introduce  three  interpretations,  of  increasing 
generality,  which  allow  us  to  associate  specific  values  to  singular 
integrals  of  the  type  under  consideration. 


a.  Principal  Value  Interpretation 


Such  integrals  were  considered,  among  others,  by  Cauchy,  who  was  the 
first  to  propose  one  interpretation,  to  which  he  gave  the  name  "Valeur 
Principale",  or  Principal  Value  ^PV),  now  commonly  denoted  by  the  letter  P. 
The  definition  is 


Principal  Value  of  Ib  = P J*  — dx 
r a Ja  x-5 


■ lim 
e*0 
£> 0 


<b  f(x,  X) 
J5+e  x-5 


dxl 


(3) 


with  a < 5 < b. 

The  conditions  for  the  existence  of  the  PV  are: 

i)  f(x,  X)  is  integrable  in  the  interval  a,h 

ii)  as  c*Q  through  positive  values,  the  limit  of  the  expression  inside 
the  curly  brackets  exists,  even  though  the  limits  of  the  individual  terms 
may  fail  to  do  so. 


The  PV  integral  of  eq.  (3)  has  been  investigated  by  Hilbert,  who  looked 
upon  it  as  the  integral  transform,  with  kernel  (x-5W,  of  f(x,X). 

Hilbert  obtained  integrals  of  this  type  while  considering  real  and  imaginary 
parts  of  Fourier  integral  transforms  of  real  functions.  Details  of  these  con- 
siderations are  given  in  Ref.  15.  Tables  of  Hilbert  transforms,  as  well  as 
various  properties  which  can  facilitate  the  extension  of  these 
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cables,  appear  in  Ch.XV  of  Ref.  14.  The  reader  should  note  ChaC  Che  PV 
integral  (3)  and  Che  formal  definicion  of  a HilberC  transform  differ  by  a 
numerical  facCor. 


Ic  should  be  nociced  thaC  having  defined  Che  PV  as  above,  we  scill  have  Che 
Cask  of  evaluacing  Che  Cwo  finice  inCegrals  inside  Che  brackets.  The 
definicion  does  not  provide  a method  for  performing  Che  integration.  But  we 
shall  see  ChaC  in  certain  cases,  where  additional  restrictions  are  imposed 
on  f(x,\),  one  is  able  Co  find  the  PV  with  sufficient  ease. 


Another  remark  that  concerns  this  definicion,  is  that  in  Che  usual 
criterion  for  Che  existence  of  l£,  namely 


lb- 

a 


lim 

rp-0 

n>0 


/r 


d« 


lim 

p*0 

p>0 


rb 


f(x.  A) 


dx 


(4) 


ri  and  p are  required  to  tend  to  zero  independently.  On  the  other  hand,  in 
eq.  (3)  we  have  restricted  the  two  limit  processes  by  setting  n ■ p ■ e. 

Hence  n,  p are  no  longer  independent  of  each  other. 

The  notion  of  PV  may  be  applied  not  only  when  the  singularity  occurs  at 
a finite  distance  from  the  origin,  but  also  when  the  integrand  FCx^)  in  the 

integral  ^ F(x,\)dx  behaves  at  + 00  in  a way  that  renders  the  integral 

divergent. 

For  this  type  of  integral  the  definition  of  PV  takes  the  form 

P £ F(x,A)  dx  = lim  f?  F(x,A)  dx  (5) 

0<R-h» 

whereas  the  ordinary  integral  of  F(x,\)  exists  if  the  integral 

lim  F(x,\)dx  (S' 

0<R-*"“  P 

0<p+«° 

where  R and  p are  independent,  converges. 

In  eq.  (5),  the  two  independent  limit  processes  of  the  ordinary 
definition  have  been  again  reduced  to  a single  limit  process. 

In  Appendix  I we  discuss  the  evaluation  of  the  PV  of  integrals  with  the 
aid  of  contour  integration  in  Che  complex  plane  for  a certain  class  of 
integrands. 

Finally  it  should  be  remarked  that  if  the  integrals  (2),  (4)  and  (6) 
exist,  then  so  do  the  appropriate  PV  integrals,  and  the  actual  values  of  the 
respective  pairs  are  identical.  The  converse  of  this  result  is  obviously 
false. 
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b.  Integral  Evaluated  Along  Deformations  of  the 
Original  Path  — The  Plemeli  Formulas 


It  is  not  possible  to  ascertain  when  and  who  was  the  first  to  propose 
path  deformation  in  order  to  deal  with  singularities  occurring  along  the 
path  of  integration.  It  will  surely  surprise  no  one  that  the  procedure 
originated  already  with  Cauchy.  After  all, Cauchy,  together  with  Riemann, 
was  the  founder  of  the  theory  of  complex  functions  and  complex  integration. 
But  the  actual  results  in  the  present  context  were  presented  much  later  by 
Plemelj,  as  described  below  in  detail. 

Since  the  method  is  widely  used  in  physics  (and  undoubtedly  in 
v engineering  as  well),  at  this  point  it  is  easier  to  give  as  reference  a 

textbook,  such  as  the  one  by  Dennery  (Ref.  1). 

. The  present  approach  consists  of  two  changes  in  our  interpretation  of 

the  original  integral.  First,  we  look  upon  the  path  of  integration  as  being 
imbedded  in  the  complex  plane,  and  second,  we  deform  the  path  so  as  to  avoid 
passing  through  the  singularity.  The  first  change  implies  that  the  real 
variable  x is  replaced  by  the  complex  variable  z * x + iy.  Concerning  the 
subsequent  deformation  of  the  path,  we  shall  make  the  further  assumption 
that  the  integrand  is  an  analytic  function  of  z.  While  this  is  not  too 
restrictive  in  applications,  it  allows  us  to  use  Cauchy's  integral  theorem 
(Ref.  2).  In  our  context,  this  theorem  implies  that  the  path  of  integration 
may  be  continuously  deformed  inside  any  region  of  the  z-plane  which  contains 
no  singularity  of  the  integrand,  without  changing  the  value  of  the 
integral.  If  such  singularities  are  encountered  during  path  deformation, 
the  new  path  must  be  indented  around  these  in  an  appropriate  fashion. 

The  procedure  is  depicted  in  Figs.  1 and  2.  The  original  path  is  shown 
in  Fig.  1,  and  Fig.  2 shows  the  two  possible  deformations  in  the  complex 
z-plane.  The  indentations  shown  there  are  smooth  curves,  but  otherwise 
quite  arbitrary.  By  Cauchy's  theorem,  for  analytic  integrands  we  may 
replace  these  arbitrary  indentations  by  semicircles  centered  at  x*fj.  This 
facilitates  the  computations  without  changing  the  actual  values  of  the 
corresponding  integrals. 


Deformed  Paths 


In  an  obvious  notation,  the  two  additional  values  of  integral  (1),  are 


proceeds  as  follows 


The  actual  evaluation  of  I..  and  I 


In  eq.  (10)  C%  denotes  the  semicircle  centered  at  x-5  and  with  radius 
e.  For  Ij,  to  exist  it  is  necessary  for  the  integrals  on  the  r.h.s.  o 
(10)  to  exist.  First  we  notice  that  the  limit  of  the  square  brackets 
represents  the  PV  of  the  integral  (1)  as  defined  by  eqs.  (2)  and  (3). 
remains  to  evaluate  the  integral  along  C£  and  then  pass  to  the  limit 


rr 


On  C*,  z * 5 + ee^,  0 £ |\J>|  <_  it.  Hence  dz  ■ ige^d^  and  we 
obtain 


j>  UzAl  iz.  i j° 

Jc  z-r  Ju  iib  y 

e® 


'C  z-F 

E * 


If  f(x,\)  is  such  that  integration  over  \|i  may  be  legitimately 
interchanged  with  the  limit  process  e-*-0,  we  obtain  after  obvious 
simplifications 

lim  J-+  dz  » -iirf(5>X) 

c*0  cc  '-5 

The  final  result  is 

lu  * / • p C 
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(12) 


(13) 


An  entirely  similar  calculation  yields 


/ 

IT 


2"C 


dz 


dx  + iirf(?,X) 


(14) 


The  above  derivation  shows  clearly  that  these  two  values  IM,  Ig  were 
obtained  by  shrinking  the  indentation,  after  the  indicated  integrations  have 
been  performed.  Thereby  I u and  Ig  have  this  limit  process  in  common  with 
the  PV  of  integral  (1).  It  is  therefore  not  surprising  that  eqs.  (13)  and 
(14)  contain,  as  part  of  the  results,  this  very  same  PV  integral. 


Another  aspect  that  should  be  highlighted  is  that  if  f(x,X)  is  a real 
function,  then  Iy  and  Ig  are  complex  functions  of  (£,X).  By  contrast, the 
PV,  if  it  exists,  is  in  this  case  a real  function  of  (£,X)« 


The  results  of  eqs.  (13)  and  (14)  were  published  first  by  I.  Plemelj  in 
1908  (Ref.  3),  in  a more  general  context.  In  his  investigations  the  path  of 
integration  was  from  the  start  a curve  on  the  z-plane,  such  that  z»£  was  on 
this  curve  and  £ was  not  necessarily  real. 


A different  representation  may  be  obtained  for  Iu  and  Ig  if  the  function 
f(z,X)  is  analytic  in  z in  the  neighborhood  of  the  real  axis.  This  is  the 
same  requirement  that  made  possible  the  choice  of  the  semicircular  paths 
C*  in  the  derivation  of  eqs.  (13)  and  (14).  In  such  a case  the  path 
is  entirely  equivalent,  by  Cauchy's  theorem,  to  a path  parallel  to  the 
x-axis  and  at  a distance  e above  it,  which  is  ultimately  let  tend  to  0. 

That  is,  the  paths  in  Fig.  3 are  equivalent. 
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Hence,  we  write 


Figure  3. 


I..  ” lim 


r®+ie 

Jj»+ie 


f(z,X) 


Now  we  change  the  variable  of  integration  as  follows: 


z * ? + iE 

Then  it  is  clear  that  -»  < C ^_®,  and  we  can  write 


I u “ lim  J_p, 


f(C+je,M 


Since  f(C,X)  is  analytic  on  the  real  axis  and  f(£  + ie,M  ^ f(?,^)  as  ^0, 
we  may  put 

. _ f®  f(x,M  , , 

IM  - lim  J*  dx  (. 

e+0  l 

£>0 


In  (18)  we  have  replaced  the  variable  of  integration  ( by  x,  only  to  be 
reminded  of  the  original  integral  (1). 

A similar  result  may  be  obtained  for  Ig : 

f®  f(x,X)  . 

- lim  J_  —Ptf  dx 


We  close  this  subsection  with  another  look  at  Itt,  Ig  as  represented  in 
eqs.  (13),  (14),  this  time  from  a different  perspective.  Since  the  function 
f(z,X)  aside  from  being  analytic,  is  completely  arbitrary,  we  can  consider 
(x-O"!  to  be  a generalized  function  (Ref.  4).  Then  the  following 
operational  interpretation  is  immediately  obtained: 


8 


(20) 


/ « ' 1 C dI  ik  - C s(x_t) 


/ ■ 'C11?;  * l*  C *■  «<«-£> 


where  (x-^)-*  and  {(x-£),  the  Dirac  {-function,  act  as  the  kernels  of  the 
indicated  integral  operators. 


We  shall  see  in  subsection  (c)  that  (20)  and  (21)  are  but  special  cases 
of  a whole  class  of  similar  results  containing  a {-function. 


If  we  remove  the  integration  signs  in  the  above  relations,  wc  obtain  the 
symbolic  formulas 

— * p " i-TT  {(x-r)  (22) 

x-£+i0  x-£  u s 


These  results  were  also  derived  in  Ref.  5,  but  in  quite  a different  way. 


In  Appendix  I we  show  how  to  evaluate  the  PV  through  the  use  of  contour 
integration  in  the  complex  plane.  First,  it  is  shown  there  that  the  simple 
requirement  of  analyticity  of  f(z,\)  in  the  neighborhood  of  the  real  axis 
has  to  be  augmented  by  more  stringent  conditions  for  the  contour  integration 
method  to  be  valid.  But  the  evaluation  of  the  PV  integral  by  this  method 
involves  essentially  the  evaluation  of  the  integrals  Iu  and  1^  (by  the  same 
method).  Thus  we  see  that  Plemelj's  formulas  (13)  and  (14)  remain  valid  for 
a larger  class  of  functions  f(z,X)  than  that  described  in  Appendix  I.  This 
remark  is  important  in  applications since  the  conditions  of  Appendix  I may 
be  too  restrictive.  Examples  to  this  effect  are  adduced  in  Appendix  II. 


The  {-Function  Interpretation 


We  return  to  the  original  integral  (1)  and  remark  that  the  integrand 
there  is  the  ratio  of  two  functions,  f(x,\)  and  x-J£.  We  may  look  upon  this 
ratio  as  being  the  solution  of  the  following  equation  in  g: 


(x-£)  g - f(x,X) 


But  in  writing  for  g the  ratio  f(x,X)/(x-£)  we  are  neglecting  a host  of 
additional  values  for  g,  which  solve  eq.  (24).  These  are  obtained  by  adding 
to  f(x,X)/(x-£)  all  the  possible  solutions  to  the  homogeneous  equation 


<x-5)  g (x,£) 


1 


At  first  sight  eq.  (25)  may  seem  to  have  only  the  trivial  solution  g0  ■ 0. 
This  is  indeed  correct  for  all  x /£.  But  for  x - £ , we  see  that  g0  can 
have  any  arbitrary  value  A(£),  where  A(£  ) may  be  a real  or  complex  function. 

To  find  the  generalized  function  interpretation  of  g0(x,£)>  we  borrow 
Lighthill's  approach  for  the  same  problem  (Ref.  6).  We  construct  from  an 
arbitrary  test  function  (x)  (Ref.  4),  another  test  function  ^(x) : 

_ 2 

^(x)  ■ 4»(x)  -ty(5)e  * / (x-O  (26) 

Notice  that  \|>i(x)  has  no  singularity  at  x^ . 

Now  we  wish  to  show  that  g0(x,O  has  the  sifting  property  of  a 
6 -function,  namely 


£ g0(x£ty(x)dx  - C(5*K5)  <27) 

where  C(£)  is  arbitrary. 

Indeed,  using  eqs.  (25)  and  (26)fwe 

£ gQ(x,5)\|)(x)dx  - (x-£)  gQ(x,£) 

- 4KO  £ gQ  (x.O 

The  last  integral  in  eq. (28)  is  seen  to  be  equal  to  some  arbitrary 
function  C(£),  since  go(x,5)  was  considered  arbitrary.  This  proves 
eq.  (27). 

If  we  choose  now  for  >|/(x)  the  unitary  function  U(x)  * 1,  for  all  x,  then 
eq.  (28)  shows  that 

£ g0(x,Odx  * C(V 

All  these  results  imply  that  we  may  write  for  go(x,£)s 

go(x,5)  - C(£)  6(x-£)  (29) 

where  6(u)  is  the  Dirac  6-function  (Ref. 6). 

The  result  of  eq.  (29)  can  be  also  proven  with  the  aid  of  Fourier 
transforms.  A discussion  of  this  approach  for  a more  general  case  is  given 
in  Section  III(c). 


find 


\|»1(x)dx  +^(^)I^  g0(x,O  eX  dx 


“X  . 

e dx 


(28) 
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Ic  is  clear  now  that 


r 


g(x)  - * c $(*-£)  (30) 

is  the  general  solution  of  eq.  (24),  for  all  x belonging  to  an  interval 
which  contains  £ in  its  interior.  If  5 is  not  inside  this  interval,  then  in 
eq.  (30)  C <5  (x-j; ) must  be  deleted.  We  also  require  f(£,X)  i 0. 

We  can  suomarize  the  preceding  discussion  with  the  statement:  division 
by  zero  is  legitimate  if  the  infinity  introduced  thereby  is  taken  explicitly 
into  account.  Having  said  this  implies  that  in  eq.  (30)  we  should  look  upon 
the  ratio  f(x,X)/(x-£)  as  being  finite,  even  at  x ■ 5 . To  accomplish  this  we 
recall  that  C is  arbitrary  and  it  will  still  be  so  if  multiplied  by  f(£,X). 
Then  we  can  write  instead  of  (30) 


g(x) 


+ C6(x-5) 


| f (x,X  ) 


(31) 


If  we  let  f(x,X)  = 1,  the  result  is  identical  with  the  one  given  in  Ref.  4. 

Now,  since  our  purpose  is  to  use  (x-£)“*  under  the  integral  sign,  we 
shall  choose  for  (x-5)”*  on  the  r.h.s.  of  (31)  a particular  determination 
that  yields  a finite  value.  In  most  cases  of  interest,  but  not  all,  the 
choice  is  the  PV  of  the  integral. 

Hence  all  possible  determinations  of  g(x)  = (x-£)“*  can  be  written 
symbolically  as 


+ C $(x-5) 
P 


(32) 


where  subscript  g and  p denote  general  and  particular  solutions, 
respectively. 


In  integral  form  we  have: 


; dx  - p ; dx  <■  Cf(5,X)  (33) 


where  we  have  taken  the  PV  for  p. 


Eq.  (32)  includes  as  particular  cases  eqs. 
includes  as  particular  cases  Iu  and  1^  of  eqs. 


(22)  and  (23),  while  eq.  (33) 
(13)  and  (14). 


It  is  of  interest  to  note  that  Jones  (Ref.  4)  defines,  in  common  with 
Lighthill,  (x-5)”l  as  that  generalized  function  which  selects  the  PV  when 
appearing  as  a kernel  of  an  integral.  Obviously,  this  is  a special  case  of 
the  general  result  in  eq.  (33). 


\ 
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III.  INTEGRALS  WITH  MULTIPLE  POLES 


Here  we  shall  consider  integrals  of  the  type 

I (5, XI  - r dx  (34) 

(x-V 

The  function  f(x,X)  has  the  same  properties  as  in  Section  II,  and  n > 1 is 
an  integer. 

a.  Principal  Value  Interpretation 

If  we  wish  to  follow  a pattern  similar  to  the  one  of  section  II  we 
encounter  a difficulty  in  attempting  to  define  the  PV  of  (34).  To  clarify 
the  situation  we  take  the  special  case  f(x,X)  = !•  With  the  usual  definition 
of  PV,  we  obtain  by  direct  calculation  for  n > 1 

r°°  dx  ■ lim  l+(*l)n  _Jo;  n odd  (35) 

~ '.-lit"1  V " ev“ 


For  odd  n the  above  result  corresponds  to  our  expectation,  since  the  area 
under  an  antisymmetric  (with  respect  to  O curve  should  be  zero.  On  the 
other  hand,  we  still  have  a divergent  integral  for  n even.  One  way  to 
remedy  this  is  to  redefine  the  PV  as  follows 


Predef 


i^iiXidx 

(x-C)n 


lim 

e>0 


[f?-€  Kiii  dx , r «*ii>  <J 
P (x-C)n  5 e (x-?)n  J 


(36) 


For  f(x,X)  = 1,  we  obtain 


j»  ® — dx — „ q for  an  integral  n > 
redefJ-  M)n 


(37) 


In  Appendix  II  we  show  that  P £ dx/(x-£)  - 0.  It  is  clear  from  (36)  that  we 


have  also 


redef 


r - 

"L®  X-Z 


(38) 


Hence  (37)  is  true  for  all  integral  n > 0. 

Whether  this  redefined  PV  is  a useful  device  or  a more  general  concept 
has  to  be  introduced,  is  still  an  open  question.  In  the  meantime, another 
possibility  exists:  we  start  with  eq.  ^34)  and  integrate  by  parts  n-1 


r 


times.  We  ass.une  that  at  x ■ + «,  the  function  f(x,X)  and  its  first  n-1 
derivatives  with  respect  to  x do  not  grow  faster  than  certain  powers  of  x, 
such  that  f^k'(x,X  J/x0"*"*1  ♦ 0,  for  0 <_  k <_  n-2. 


Then  we  find 


Z 


«»A)  dx 


(n-1 )! 


f(°~U(».X) 


dx 


(39) 


Using  this  relation,  when  the  conditions  of  its  validity  are  satisfied,  we 
define  the  PV  of  (34)  as 


9 


f(x,X) 

(x-on 


dx  i 


(n-1)! 


-(n-1).  , \ 

p£  f JgjAJ  dx 


T^TJT 


lim 

e>0 


[* 


- ( O” 1 ) / X X 

f (x,X ) 


dx  + 


r 


X< 


f(n~U(x,X) 


dx 


] 


(40) 


We  have  used  script  P to  denote  this  more  general  type  of  PV. 


We  note  that  with  this  definition,  if  f(x,X)  5 1 for  -®  <_  x <^«,  we 
obtain  the  same  result  as  in  eq.  (37). 

In  eq.  (39)  we  have  assumed  that  all  the  derivatives  concerned  exist. 

If  we  interpret  the  function  f(x,X)  to  be  a generalized  function,  such  as 
defined  in  Refs.  (4)  and  (5),  these  derivatives  always  exist  and  are 
generalized  functions  themselves.  When  we  replace  f(x,X)  in  (39)  and  (40) 
by  [H(x-a)  -H(x-b)]  f(x,X)  we  obtain  the  result  given  by  Jones  in  § 4.4  of 
Ref.  4,  since  during  the  integration  by  parts  procedure,  the  integrated 
terms  no  longer  vanish.  In  the  same  reference  one  finds  the  connection,  for 
this  case,  with  the  "Finite  Part"  concept  introduced  by  Hadamard  (Ref.  11). 
Jones  (Ref.  4)  seems  to  have  been  anticipated  by  Fox  (Ref.  12),  who  used 
Hadamard 's  ideas  to  generalize  the  PV  of  an  integral  of  type  (34),  along  a 
finite  interval.  Fox  obtains  basically  the  same  result  as  eq.  (40)  for  the 
infinite  interval,  as  well  as  the  generalization  of  Plemelj  formulas  given 
below.  It  appears  that  the  derivations  of  Fox,  based  on  standard  analytic 
arguments,  apply  to  a more  restricted  class  of  functions  than  the 
interpretation  in  terms  of  generalized  functions  given  by  Jones. 

b.  Plemelj* s Generalized  Formulas 


By  analogy  with  the  treatment  for  n ■ 1,  we  define  two  integrals 
along  deformed  paths: 


f(x,X) 

on 


dx 


(41) 
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4BNBK 


5 r 


(42) 


The  paths  in  (41)  and  (42)  are  those  shown  in  Fig.  2.  The  difficulty 
encountered  in  (a)  above  is  present  here  too  when  attempting  to  evaluate 
(41)  and  (42)  according  to  the  prescription  given  in  eq.  (10).  This  becomes 
clear  when  we  write  down  the  equation  corresponding  to  eq.  (11): 


r ) 

K (^>0 


dz  ■ 


_1_  r° 

n-1  -it 


e-i(n-l>j)f(^en[)  A)<4J> 


(43) 


It  is  seen  immediately  that  when  e*0  the  r.h.s.  of  eq.  (43)  diverges. 
Moveover,  the  contribution  to  Iy  from  the  path  along  the  real  axis  also 
diverges,  as  we  have  seen  in(a). 


Therefore  we  proceed  here  as  in  (39)  and  first  write 


f liiAl 

J%(*-{-)n 


dz  ■ 


dz 


(44) 


On  using  the  prescription  of  eq.  (10)  for  the  r.h.s.  of  eq.  (44),  we 
obtain 


i - r — dz 


PZ 


tUjXl 

(x-^)n 


dx 


(n-1)! 


In  a similar  fashion  the  analogue  of  eq.  (14)  is 


(n-1), 


y / - ■ PC  f * * *4?^ 

* (z'5)  (x^} 


(45) 


(46) 


"V 

We  note  that  (45)  and  (46)  reduce  to  Plemelj's  formulas  (13)  and  (14) 
when  n ■ 1.  Hence  we  call  (45)  and  (46)  Plemelj's  Generalized  Formulas. 


The  following  results  correspond  to  the  results  in  eqs.  (18),  (19), 
(20),  (21),  (22)  and  (23): 


I - 
U 


lim 

<~0 

E>° 


1 

(n-1)! 


x-^+ie 


dx 


f(n~0(x.\) 

x-^-ie 


dx 


(47) 


(48) 


"/V  (*^ ) 
! ^ 


-v*  (*-s 
i 


PSZ 

dx 

fcK-l)0-1 

■£  dx5(n“U  (x-5) 

(x-5  )“ 

(n-1) ! 

PC 

dx 

. in  (-I)0-1 

■C  d*s(n_1)  (x-5) 

(x-S)n 

(n-1 ) ! 

p- 


(x-^+iO)°  w (x-0° 


iirC-l)0"1  , (n-1)  /„  ^ 

(n-iri'  5 (jK) 


(49) 


(50) 


(51) 


— * — p — — * a (-i)”;1  s<-‘>ch.) 

(x-C+iO)n  (x*On  (n-U! 


(52) 


c.  The  6 -Function  Interpretation 

In  this  subsection  we  are  concerned  with  the  general  solution  to  the 
equation 


(x-£)ng(x)  ■ 1 

For  n ■ 1,  we  have  found 

«l<*>  ’^5  | 


(53) 


(54) 


where  the  subscript  p denotes  a particular  interpretation  of  (x-C)_1, 
which  in  eq.  (32)  we  took  to  be  P,  the  Principal  Value. 

« 

We  proceed  by  induction.  Having  solved  for  g^(x)  from  the  equation 

(x-C)8l(x)  - 1 (55) 

we  solve  now  for  g^(x)  from 

(x-S )2g2(x)  - 1 (56) 

But  we  can  write  (56),  on  using  (55),  as 

(x-?)g2(x)  « gj(x)  (57) 

Use  of  (54)  on  both  sides  of  (57)  yields, 


(x-cr 


+ A. 
P 1 


fi(x-C) 

X-C 


♦ A2  6(x-5) 


(58) 
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To  bring  (58)  to  the  final  form, we  prove  the  following  result  concerning  the 
6 -function  and  its  derivatives: 


6(r)(x)  - -rf  (r+1)(x) 


(59) 


for  all  integral  0 <_  r. 


Multiplying  (59)  by  an  admissible  function  and  integrating  by  parts  its 
r.h.s.,  it  is  easily  proved  that  (59)  is  indeed  valid. 

From  (59)  we  obtain 


g(r)(«-g)  . _ 6(r+l) 
xs  " 


(x<) 


(60) 


and  in  particular  for  r ■ 0,  we  have 


— - - « \x<) 


(61) 


Hence  we  can  write  for  (58) 


8o(x>  5 


1 


- — 7 B + c6(x-5)+c/(x-5) 

i.  p o 1 


u<r  18  (x-n 

where  Cq»  are  arbitrary  constants. 

It  is  clear  now  that  proceeding  in  this  way  and  using  eq.  (60) 
repeatedly  we  find  for  arbitrary  integral  n > 0 

— - — I + “j1  C <5(r)  (x-5) 

(x-5)n  8 (x-f;)n  1 p £ r 


r-0 


(62) 


(63) 


Without  specifying  the  particular  value  p,  this  result  is  given  by  Jones 
in  eq.  I 6.2  of  Ref.  4. 

In  eq.  (63)  we  choose  for  p the  PV  as  defined  in  eq.  (40).  Hence  we  can 
finally  write  the  symbolic  formula 


i-r  | * “j1  c «(r)<«-5) 

.Mn  I 8 (x<)n  Z r 


(x-5)‘ 


r-0 


(64) 


Using  this  result  we  can  write  the  analogue  of  eq.  (33) 


£ dx  . PZ  . "Ji  (r) 

(X-C)n  (x-On  rt0  Kr  f 


(65) 


16 


The  general  solution  of  eq.  (53)  presented  in  eq.  (64)  ’ ■?.  'Obtained 

also  by  applying  a Fourier  transform  to  eq.  (53),  solving  the  resulting 
linear  differential  equation  in  G(k)  (the  transform  of  g(x)),  and  then 
performing  the  inverse  Fourier  transform  which  yields  g(x).  The  arbitrary 
constants  Cr  make  their  appearance  as  the  constants  of  integration  of  the 
differential  equation.  Since  the  differential  equation  is  inhomogeneous 
(i.e.,  it  has  a non  zero  term  on  the  r.h.s.)  its  solution  consists  of  the 
sum  of  a particular  solution  and  the  general  solution  to  the  homogeneous 
equation.  This  is  precisely  the  result  presented  in  eq.  (63). 


IV.  HALF-RANGE  INTEGRALS 
In  this  section  we  consider  integrals  of  the  type 

I(X,5)  -T0  dx  (66) 


where  £ > 0.  This  is  representative  of  the  half-range  integrals  along  the 
negative  real  axis  as  well,  since  a change  of  variable  will  bring  the 
integral  to  the  form  (66),  provided,  of  course,  that  f(x,X)  has  the  required 
properties.  Furthermore,  integrals  of  the  type 

I(X,5,a)  - dx  (67) 

where  5 > a >_  0,  are  expressible  in  the  form  (66)  by  a translation 

I(X,5,a)  - I(X  ,5 -a)  - ro  vgl")  dx  (68) 

where,  by  the  previous  conditions,  we  have  EJ  - a > 0. 

It  is  not  difficult  to  show  that  Plemelj  formulas  (13),  (14)  are 
applicable  here  as  well,  such  that 

Iu  = J dz  - P /J  ~^X—  dx  - inf(C,X)  (69) 

Y ; dz  - P ro  dx  + i*f(5,X)  (70) 

The  same  is  true  of  the  6-function  interpretation  as  discussed  in  11(c). 

Therefore  it  remains  to  evaluate  the  PV  of  integral  (66),  by  any 
method.  First,  we  note  the  following 


5>o 


(71) 


Hence,  unlike  the  full-range  integral,  the  PV  of  (66),  for  f(x,X)  = l, 
diverges.  This  result  immediately  precludes  the  evaluation  of  the  PV  bv 
subtraction  of  the  singular  part  as  shown  in  (1-17)  for  the  full-range 
integrals.  Also,  writing  f(x,X)  as  a sum  of  an  odd  and  an  even  function 
does  not  yield  any  advantage  in  the  present  case.  Nevertheless,  several 
possibilities  do  exist  for  the  PV  evaluation  of  (66),  depending  on  the 
nature  of  the  function  f(x,X). 

A finite  particular  value  to  be  associated  with  the  integral 
Jo  dx/(x-C)  is  offered  in  Appendix  VI,  in  the  spirit  of  the  Finite  Part 
of  Hadamard  (Ref.  11)  and  the  6-function  interpretation  of  Section  11(c). 

a.  PV  of  Full-Range  Integral  is  Known 

If  we  assune  that  f(x,X)  is  defined  over  the  full-range  |x|  < ®,  and 
that  the  PV  can  be  explicitly  evaluated,  i.e., 

p Si  dx  M F(C,X)  f 72) 

then,  if  C > 0,  we  can  write 

P Jq  dx  - F(C,X)  ♦ J*  dx  (73) 

The  integral  on  the  r.h.s.  of  (73)  is  an  ordinary  integral  now  (the 
integrand  has  no  singularities)  and  moreover  it  is  a convergent  integral, 
since  F(£,X)  is  finite  by  assumption.  Thus,  formula  (73)  allows  us  to 
dispense  with  the  t limit  in  the  PV  definition.  This  approach  may  be 
practical,  since  the  form  of  the  integral  on  the  r.h.s.  of  (73)  shows  it  to 
be  the  Stieltjes  transform  of  f(-x,X).  Tables  of  such  integrals  are  given 
in  Ref.  ( 14). 

b.  Evaluation  by  Contour  Integration 

Let  us  assume  that  the  appropriate  requirements  are  satisfied  by  f(x,X^ 
when  extended  to  the  sector  of  the  complex  plane  given  by  |argz|  < a, 

0 < a < ir/2. 

To  evaluate  the  PV  of  ^66)  we  basically  evaluate  Iw  and  I in  a fashion 
entirely  analogous  to  the  one  described  in  Appendix  I.  But  here  the  results 
are  not  the  same,  as  it  will  become  clear  from  the  typical  case  discussed 
below. 

First,  we  evaluate  IM,  by  choosing  the  contours  in  Figs.  4 and  5,  each 
contour  corresponding  to  the  specific  behavior  of  f(z,X). 


Fig.  4 


Fig.  5 


Thus  for  X > Xo  the  contour  of  Fig.  4 may  be  the  appropriate  one.  We 
assume  that  the  integral  along  ^vanishes  in  the  limit  R •+■  ®,  and 
the  final  result  is 


tr0  ^ - mcc.x>  * j - 


y-le 


***■  0<ar£2<a  Res  W1'  X<X0 


(74) 


Similarly,  for  X < Xot  the  contour  of  Fig.  5 may  be  the  appropriate  one, 
and  with  the  integral  along  ^ vanishing  with  R ■+•  oo,  we  find 


f (x,X) 
*-£ 


dx 


iirf (C  ,X)  + f(g  dy 


y~Ce 


27T1-a<ar|z<0 


Res 


z-  £ ’ 


X<Xr 


(75) 


We  observe  immediately  that  the  integrals  on  the  r.h.s.  of  (74)  and  (75)  are 
ordinary  integrals,  devoid  of  singularities  along  the  path  of  integration 
(provided  a is  so  chosen  that  f(z,X)  has  no  singularities  along  the  ravs 
z - ei.^;  0 < v). 


The  evaluation  of  I*  proceeds  in  exactly  the  same  way  and  we  therefore 
omit  it. 


c.  Evaluation  by  Laplace  Transform  Identity 
We  make  use  in  this  section  of  the  identity 

jT  e au  du  * — ; Re  a >0  (76^ 

J 0 a 
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This  is  none  ocher  than  Che  Laplace  transform  of  the  Heaviside  function  H(u). 

Now  we  can  write  from  the  definition  of  PV  and  of  an  integral  with  an 
infinite  upper  limit: 


| pit*  HP- 

Hence,  on  using  (76)  we  obtain 

z5"  ax  - -/En-cdxf(«,Mn  ;<5-x)u  du 

0 x~^  u u 

• -J^  du  e~^u  Jg^dx  f(x,M  eUX 


dx 


(77) 


(78) 


where  we  have  assumed  that  the  interchange  of  integrations  is  permissible. 
Similarly  we  find 

J?4£  dx  - ^ do  e*  f^dx  £<x,X)e'UX  (79) 

On  adding  (78)  and  (79)  and  taking  the  limits  indicated  in  (77),  we  find 

P d,  • do  le"5  dx  - o'1*  l\  dx  f(x,X)  .“|  (80) 

The  integrals  appearing  in  (80)  are  ordinary  integrals  (the  integrands  have 
no  singularities  along  the  paths  of  integration)  which  converge  if  f(x,M 
has  the  appropriate  behavior  as  x ®. 


Ue  remark  that 
fail.  An  example 
will  be  discussed  in  Appendix  IV. 


the  present  method  may  work  when  methods  (a)  and  (b)  will 
of  such  behavior  arises  for  f(x,M  " e“^x,  ^ > 0,  which 


Incidentally,  method  (c)  may  be  advantageous  in  evaluating  the  PV  of 
full-range  integrals,  although  we  have  not  used  it  for  this  purpose  in  the 
present  report. 

We  consider  now  the  case  of  higher  order  poles.  First,  it  is  easily 
seen  from  the  definition  that 


dx 


(x-5)1 


(-1) 


n-1 


(n-l)C 


n-1 


n even 
n odd 


(81) 


where  n > 1. 
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Hence,  as  in  Che  case  of  Che  full-range  inCegral,  we  musC  generalize  Che 
definicion -of  PV  in  Che  same  way  ChaC  was  done  in  SecCion  III. 

The  definicion  is 


f (xA) 


0 (x-Ol 


__1 

(n-1) ! 


P r f(n~1)(^> 

r J n 


We  can  now  make  use  of  all  Che  meChods  known,  such  as  (a),  (b)  and  (c) 
above,  Co  evaluaCe  Che  r.h.s.  of  (82). 

The  Plemelj  formulas  may  be  obCained  now  by  replacing  Che  PV  inCegral  on 
Che  r.h.s.  of  (82)  wich  Che  appropriate  expressions,  as  is  done  in  (45)  and 
(46). 


V.  INTEGRALS  WITH  FINITE  LIMITS 

The  inCegrals  we  discuss  here  are  of  Che  type 

1(5  ,A)  ■ -^X-}  dx  (83) 

where  a < 5 < b. 

The  PV  is  defined  as  in  eq.  (3).  This  PV  can  be  considered  as  Che 
finiCe  Hilbert  transform  (up  Co  a numerical  factor).  A few  references  to 
this  type  of  PV  inCegral  may  be  found  in  Ch.  XV  of  Ref.  14.  The  acCual 
evaluation  of  the  PV  is  hampered  by  our  inability  to  find  the  anti-derivative 
of  Che  inCegrand.  Therefore,  we  again  musC  resort  to  other  methods  in 
order  to  avoid  Che  e-limit  process. 

But  if  the  PV  of  (83)  exists,  then  the  Plemelj  integrals  IM  and  1^  also 
exist:  * 

V / dz  » p !ha  dx  - iTTf(5,A)  (84) 

I^i  j dz  - P /b  dx  ♦ inf (5 ,A)  (85) 

The  6-function  interpretation  is  clearly  still  valid  in  the  present  case 
also,  and  the  pertinent  formulas  carry  over,  with  minor  changes,  from  11(c). 

Using  definition  (3),  we  find  immediately 

P ^ ^5  “ lo8  * a<5<b  (86) 

We  list  now  several  methods  of  evaluating  the  PV  of  (83). 
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The  inCegrel  appearing  in  Che  r.h.s.  of  (88)  is  again  an  ordinary  convergent 
integral.  It  is  not  difficult  to  verify  that  the  expression  obtained  when 
multiplying  the  denominator  by  its  complex  conjugate  has  no  zeros  on  the 
real  axis. 


c.  Evaluation  by  Laplace  Transform  Identity 
Use  of  (76)  in  the  present  case  yields 

P dx  - £ du  I*"5  Ji  dx  f(x,X - e--5  S\  dxfCx.X)."*  | (89) 


Ue  notice  that  the  inner  integrals  are  independent  of  the  behavior  of  f(x,X) 
at  infinity,  but  the  convergence  of  the  integral  over  u may  depend  on 
f(x,X).  For  the  convergent  case,  this  formulation  allows  us  again  to 
dispense  with  the  e-limit  process. 


d.  Evaluation  from  PV  of  Full-Range  Integral 


This  method  will  apply  when  the  function  f(x,X)  can  be  extended  along 
the  entire  real  axis  and  the  PV  of  the  full-range  integral  exists.  Then  we 
can  write 


p UlAi  dl  - p £ dI  - £ 1<2^>  d,  - £ dx  (90) 


Since  we  have  assumed  a < 5 < b,  the  last  two  integrals  are  ordinary 
integrals  which  moreover  converge,  since  the  PV  integral  over  the  whole  real 
axis  is  assumed  to  exist. 


To  treat  poles  of  higher  order,  we  follow  the  procedure  mentioned  in 
Section  111(a).  The  resulting  formula,  after  repeated  integration  by  parts, 
is  more  complicated  since  we  have  to  keep  track  of  the  contributions  from 
the  endpoints  of  the  interval  of  integration. 


The  result  is  therefore  the  generalized  form 

r.(n-l) 

a 


^>;b  ilsAi  d,  . 1_  p /*>  £' 

* (x-5)“  <n'UI  ' 


1 


(n-1) ! V (n“2“k)! 


x-S 
(k) 


(x.X) 


dx 


(x.X) 


k-0 


(x-5) 


n-l-k 


b 

a 


(91) 


The  Plemelj  integrals  Iy,  I#  for  this  case  are  now  obtained  in  the  same  way 
as  in  (45)  and  (46),  by  replacing  the  PV  integral  on  the  r.h.s.  of  (91)  with 
the  appropriate  expressions. 
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VI.  CONCLUSION 


The  present  report  has  shown,  in  detail,  hew  to  interpret  integrals  with 
poles  located  on  the  path  of  integration.  In  particular,  through  the  use  of 
concepts  from  the  theory  of  generalized  functions,  it  was  established  that 
such  integrals  possess  an  infinite  set  of  values. 

The  choice  of  a value,  to  be  considered  the  "correct"  value,  depends 
entirely  on  the  context  from  within  which  the  singular  integral  arose.  The 
ability  of  making  such  a choice  depends  in  turn  on  knowing  the  set  of  values 
associated  with  a particular  integral.  It  was  shown  that  knowledge  of  the 
Principal  Value,  when  it  exists,  is  equivalent  to  knowledge  of  the  entire 
set  of  values.  Thus,  several  practical  methods,  which  replace  the  limit 
process  appearing  in  the  definition  of  the  Principal  Value,  have  been 
presented  for  the  various  cases  under  consideration  (full  and  half-range, 
and  finite  intervals). 

Several  examples  were  worked  out  in  detail,  to  illustrate  the  theory  and 
to  exhibit  the  usefulness  as  well  as  the  limits  of  applicability  of  the 
various  methods  of  actual  evaluation. 

A feature,  of  both  practical  and  mathematical  interest,  is  that  the 
uniqueness  of  the  Principal  Value  yields  integral  identities.  Such 
identities  enable  us  to  substitute  a rapidly  convergent  integral  for  a 
slowly  convergent  one.  The  advantages,  numerical  and  otherwise,  of  such  a 
substitution  are  well  known. 

The  main  purpose  of  the  work  reported  here  was  to  gather  all  the  known 
results  on  integrals  with  singularities,  fill  the  gaps  where  necessary,  and 
present  the  whole  subject  in  a coherent  and,  hopefully,  a usable  form.  As 
the  work  advanced,  it  became  clear  that  certain  omissions  will  be  necessary 
in  order  to  keep  the  report  at  a reasonable  length.  Several  important 
omissions  that  deserve  a separate  treatment  are:  (i)  integrals  with  branch 
point  singularities,  (ii)  integral  with  a singularity  (pole  or  branch  point) 
occuring  at  one  endpoint  of  the  interval  of  integration,  and  (iii)  singular 
integrals  in  spaces  of  higher  dimension.  The  author  hopes  to  return  to 
these  subjects  in  future  reports. 

A final  remark  concerns  the  theory  of  singular  integral  equations  where 
PV  integrals  of  the  type  discussed  here  make  their  appearance.  The  function 
f(x,X)  in  this  case  is  not  given  and  must  be  found  as  a solution  of  the 
integral  equation.  Fundamental  work  in  this  field  is  due  to  S.  G.  Mikhlin, 
whose  contributions  led  later  on  to  the  theory  of  pseudo-differential 
operators,  a subject  of  current  interest.  Ref.  16  contains  an  appraisal  of 
this  work  and  extensive  references. 


APPENDIX  I 


EVALUATION  OF  THE  PRINCIPAL  VALUE 

The  PV  of  an  integral  over  a finite  interval  a,b  was  defined  in  eq. 
(3).  When  we  take  the  additional  limits  a •+■  b •+■  +«  independently,  the 
integral  may  or  may  not  converge.  In  both  cases  we  shall  denote  by 
Principal  Value  of  a singular  integral  the  expression  obtained  from  the  two 
limit  processes 


f(x,  X) 


J-oo  x-j; 


lim  lim 
R-*cd  g+0 
£>0 


frC“E  f(x,  X) 

in  v~r 


dx  + ]* 

JC+E  X-C 


(1-1) 


taken  in  the  indicated  order.  Thus  (1-1)  may  exist  when,  by  the  usual 
definition,  the  integral  diverges. 

The  evaluation  of  the  PV  requires  knowledge  of  the  primitive  of 
f (x,  X)/(x-{;) . But  it  is  well-known  that  no  general  method  exists  for 
determining  the  anti-derivative  of  an  arbitrary  function.  Hence  the 
analytic  evaluation  of  the  PV  is  no  easy  task.  Fortunately,  for  many  cases 
of  interest  in  applications,  the  function  f(x,X),  when  extended  to  the 
complex  z-plane,  has  properties  that  enable  us  to  obtain  the  PV  as  a 
by-product  of  contour  integration. 

The  method  is  best  explained  when  we  refer  to  Fig.  7 below.  First,  we 
assume  that  f(z,X)  is  analytic  in  the  upper  half  of  the  z-plane.  Then  we 

take  the  integral  around  the  closed 
path  shown  in  Fig.  7.  If  f(z,X)  has 
^ singularities  inside  this  contour, 

/ \ ^ such  as  a number  of  simple  poles  zy  f 

/ \p  then  Cauchy's  theorem  yields 


Figure  7 


,f(Xj_X).  dx  + jR  dx  + f — dz 

ar  x-£  je+e  x-r  J_+  z-r 


'C+E  X-C 


+ Jc+  %!^X)  dz  ■ 2iri  l 


f(z,X) 

z-C 


(1-2) 
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This  equation  is  valid  as  it  stands,  the  only  requirement  being  the 
analyticity  of  f(z,X),  as  imposed  above.  For  this  result  to  be  useful  in 
the  evaluation  of  the  PV,  the  function  f(z,A)  must  satisfy  two  additional 
requirements:  (i)  f(£,X)  is  finite;  (ii)  the  integral  along  the  semicircle 
C+  must  vanish  in  the  limit  E ■*  •. 


Note  that  we  have  already  demanded  that  (i)  be  satisfied  in  our  initial 
discussion  of  the  Principal  Value. 


To  see  what  requirement  (ii)  entails,  we  write  for  the  integral  along  C+ 


f „ dz  - « f»  .*«d« 

V Jo  IUl9-5 


. i j-ir  d9 

» 1 "I  •‘19 


(1-3) 


Hence 


f „ dz  I 

J z-r 


1 -1 


.-l9 


f -10,  W 


_1  „ lf(8.i9,>)ld 


(1-4) 


It  is  clear  from  (1-4)  that  (ii)  implies  then  the  requirement 

fii-a)  | f ( z , x ) J •+•  0 as  | z|  + ao,  uniformly  in  arg  z,  for  0 arg  z » 0 ^ n. 

This  condition  is  the  modified  form  of  a requirement  placed  on  f(z,\),  in  a 
slightly  different  context,  by  Whittaker  and  Watson  (Ref.  7 — § 6.22).  A 
particular  case  of  (ii-a)  is  a function  of  the  form  f(x,^)  ■ eLAxg(x). 
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Then  if  g(z)  is  analytic  and  X > 0,  the  condition  (ii-a)  on  f translates 
into:  |g(z)|  * constant  as  |z|  ”*■  °°,  uniformly  in  arz  z,  for  0 ^ arg  z < it. 
This  is  a modified  form  of  Jordan's  lemma  (Ref.  7 - * 6.2221. 


We  finally  obtain,  when  taking  the  limits  (indicated  in  eq.  (I-D)  in 
eq.  (1-2) 


P C dx  - iTTf(E,X>  - z*i  Z Res 


f(z,X) 


z»z. 


(1-5) 


The  second  term  on  the  l.h.s.  of  (1-5)  is  the  contribution  of  the 
integral  along  c£  (when  t •*  0),  while  the  sum  is  extended  over  the 
residues  of  all  poles  in  the  upper  half  plane. 


Since  the  PV  integral,  if  it  exits,  is  unique,  its  value  should  not 
depend  on  the  manner  of  its  evaluation.  Hence,  if  the  semicircle  of  radius 
£ is  taken  below  the  real  axis,  we  should  obtain,  all  other  things  being 
equal,  the  same  value  as  in  eq.  (5). 


From  Fig.  8 we  can  write  down 


The  appearance  of  f(£  ,X  ) on  the  r.h.s.  of  (I~6)  is  due  to  the  fact  that 
the  closed  contour  includes  now  the  pole  at  z * £ . On  the  other  hand  we 
have 

f dz  - i J*  dd)  fQJ+ce^A) 

C~  ° 

e 

Therefore,  when  E + 0 we  obtain 

lim  / * iiTf(![  ,X) 

e+0  C“ 

€>0  e 

Collecting  all  these  results,  we  find 

P £ dx  ♦ iirf<5,X)  - 2ni  jf(^,X)  ♦ E 

If  the  function  f(z,X)  is  not  analytic  in  the  upper 
not  satisfy  fii-a)  there,  we  may  close  the  path  of  integration  with  a 
semicircle  in  the  lower  half  plane.  We  show  how  the  form  of  the  result  is 
affected  by  this  change. 

Fig.  9 deoicts  the  closed  contour  in  this  case. 

Here  requirement  Cii-a)  is  changed 
into:  (ii-b)  |f(z,X)|  ♦ 0 as 

|z|  ■*■  •,  uniformly  in  arg  z,  for 
IT  < arg  z * 9 ^ 21T.  For  f(x,X)  * 
e"^x  g(x),  with  X > 0 and  g(z^ 
analytic,  condition  (ii-b'i  on  f 
translates  into:  |g(z)|  * constant 
as  |z|  ■*■  *,  uniformly  in  arg  z,  for 
IT  < arg  z < 2tt.  Again,  this  is  a 
modified  form  of  Jordan's  lemma. 

Fig.  9 


Since  the  direction  of  integration  here  is  opposite  to  that  of  the 
previous  case,  being  in  fact  the  opposite  to  the  conventional  positive 
orientation  of  the  contour  used  in  Cauchy's  residue  theorem,  we  obtain 


E**kJ 


(1-7) 


half  plane,  or  does 
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fOtjXl 

x-C 


dx  * 


d z + 


2iri 


f (5  ,X  ) + Z Res 

j 


f(z.X ) 

2-5 


Z“Z  . 

J J 

(1-8) 


Using  (ii-b)  end  teking  the  limits  e ♦ 0 end  R ♦ • yields  the  result  for 
this  cese 


p C f^|X)  d*  - iwf(5,X)  - - 2iri 


f(S,X)  ♦ £ Res 

j 


z*z . 


(1-9) 


An  identicel  result  may  be  errived  et,  if  insteed  of  the  upper 
semicircle  c£,  we  choose  the  lower  semicircle  C£.  In  this  cese 
the  closed  contour  does  not  include  the  pole  et  z ■ £,  but  then  C£ 
gives,  for  e *►  0,  the  correct  contribution. 

It  is  elso  important  to  remark  thet  the  two  PV's  appearing  in  eqs.  (1-5) 
and  (1-9),  should  not  be  identified,  since  those  two  belong  to  different 
functions.  One  has  the  requisite  properties  in  the  upper  half-plane  while 
the  other  possesses  them  in  the  lower  half  of  the  z-plane. 

We  notice  that  the  l.h.s.  of  eqs.  (1-5)  and  (1-9)  represent  precisely 
Plemelj's  integrals  Iu  and  Ig  shown  in  eqs.  (13)  and  (14).  Therefore  we 
can  sunmarize  all  of  the  above  results  as: 


p C dx  " 


iuf(5,X)  ♦ 2iri  Z Res  f(V° 
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z*z. 


for  f satisfying  ii-a) 


- iirf(5,X)  - 2*i  £ Res 
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for  f satisfying  ii-b) 


2iri  Z Res  f(z.X) 

k z-C  z-z. 


for  f satisfying  ii-a) 

- »i  . t »..  ] 

for  f satisfying  ii-b) 


(1-10) 


(1-11) 
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2ni  f(C,X)  + Z Res  f(z.X) 
k z-C 


z-z. 


for  f satisfying  ii-a) 


l - 2iri  Z Res  Oli-P 

i 1-5 


z*z . 
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for  f satisfying  (ii-b) 


(1-12) 


It  is  important  to  observe  that  Plemelj's  integrals  Iu  and  I|,  which 
originally  were  defined  along  indented  paths,  with  the  indentations 
ultimately  shrunk  to  zero,  can  be  evaluated  by  contour  integrations 
(whenever  the  appropriate  conditions  are  satisfied)  as  if  the  indentations 
were  fixed.  That  this  is  true  can  be  immediately  seen  from  the  r.h.s.  of 
eqs.  (1-11)  and  (1-12).  No  reference  is  made  there  to  the  PV's  of  the 
integrals  (although  these  PV's  enter  the  result  implicitly).  Hence  IM  and 
1^  can  be  interpreted  also  as  integrals  along  indented* but  fixed,  paths 
and  their  evaluation  performed  by  contour  integration.  Such  integrals, 
which  avoid  the  singularities  in  their  integrands,  are  familiar  from  the 
theory  of  the  inverse  Laplace  transform  (where,  however* t he  paths  separate 
the  z-plane  into  left  and  right  halves). 

A further  remark  concerns  the  frequently  met  case  when  f(x,  X)  is  a real 
function.  It  may  happen  then  that  the  extension  of  f(x,X^  to  the  complex 
plane,  namely  f(z,X),  does  not  satisfy  either  one  of  the  requirements  (ii-a) 
or  (ii-b).  Under  these  circumstances  it  would  appear  that  the  PV  cannot  be 
evaluated  by  contour  integration.  Nevertheless,  this  difficulty  may  be 
circumvented  if  it  is  possible  to  find  a complex  function  F(z,X)  which 
satisfies  at  least  one  of  these  requirements,  and  such  that  f(x,  X)  is  equal 
to  the  real  or  imaginary  part  of  F(z,X)  in  the  limit  y -*•  0.  Two  well-known 
examples  are  fi(x,X)  * sin  Xx  and  f2(x,X)  ■ Jn(  Xx) . The  appropriate 
functions  are:  F^(z,X)  * exp  (iXz)  with  fj(x,  X)  * Im  {Fi(x,X)},  and 
F2(x,X)  ■ H^D(Xz)  with  f2(x,  X)  * Re  {^(x,  X)  }.  Use  of  F^  and 
F2  in  the  contour  integration  will  ultimately  yield  the  desired  results 
for  fj  and  f2« 

When  contour  integration  is  not  convenient  or  impossible,  we  can  go  back 
to  the  original  definition  of  the  PV.  Such  cases  may  occur  when  f(z,  X)  has 
branch  points  in  the  z-plane  or  when  it  is  not  analytic.  An  example  for  the 
latter  case  is  f(z,X)  * exp(-X|z|),  X > 0. 

A simple  procedure  is  useful  in  obtaining  the  PV  exactly  or,  at  the  very 
least,  an  alternative  expression  for  it  that  does  not  involve  the  limit 
process  t ■*  0.  This  is  based  on  the  fact  that  any  function  g(x)  may  be 
represented  as  a sum  of  two  functions,  one  odd  and  the  other  even: 

g(x)  * ■|[’g(x)  - g(-x)^  + -|Cg(x)  + g(-x7)  (1-13) 
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To  use  this,  we  change  the  variable  of  integration  in  the  integral  to 
obtain 


C d*  ■ £ d- 


Cl-IA'* 


Identifying  g(x)  with  f(x+£,X)  we  find 

r ««•»  d, . r «»*s.x)  - iiaiSiii  dx . 

■L«  x-E  ■*-<»  2x 


+ £ f(«*C»X)  ♦ f(-x+C,X)  dx 


(1-15) 


The  first  integral  on  the  r.h.s.  of  (1-15)  has  an  even  integrand.  Hence, 
if  this  integral  exists,  it  dpes  not  vanish.  Moreover,  the  integrand  tends 
to  f (5)  when  x 0,  and  if  f (£)  is  finite,  there  is  no  singularity  at 
the  origin.  Therefore, with  appropriate  behavior  of  f(x,X)  at  x • + «,  the 
first  integral  will  exist  in  the  conventional  sense.  On  the  other  hand, the 
integrand  in  the  second  integral  has  a singularity  at  x ■ 0,  but  because  the 
integrand  is  odd,  the  definition  of  the  PV  yields  zero  for  this  integral. 

Finally,  we  may  write 

p f H«,X)  d . r-  dx 

x-E  2x 

Apart  from  a numerical  factor,  eq.  (1-16)  was  obtained  by  Hilbert  (Ch. 

V of  Ref.  15). 


(1-16) 


Thus  we  see  that  the  PV  is  given  bv  a conventional  integral,  which  mav 
still  be  difficult  to  evaluate,  but  at  least  the  limit  process  e -*■  0 has 
been  therebv  eliminated.  The  same  is  true  for  Plemelj’s  formulas  '13)  and 
'14). 

A method  somewhat  similar  to  the  above  is  to  subtract  from  (1-0  the 
singular  integral  P £ f(C*X)/(x-£)(dx),  whose  value  is  obviously  zero. 

Then  we  obtain  on  the  r.h.s.  a conventional  integral 

, r t, . r --Wtq)  dt  n.m 

x-fc,  *•  x-e 

provided  f’(£,X)  is  finite  for  the  values  of  £ and  X considered. 
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appendix  II 


LIST  OF  FULL-RANGE  SINGULAR  INTEGRALS 

In  this  Appendix  we  present  a list  o£  the  values  that  can  be  associated 
with  the  singular  integrals  of  several  functions.  This  list  is  constructed 
by  showing,  in  some  detail,  how  to  apply  the  theory  developed  in  the  text  to 
the  various  functions  under  consideration.  Here  the  particular  values  used 
are  all  PV's  and  the  additive  term  of  the  ^-function  interpretation  are 
understood,  but  not  shown. 


1.  f(x,X)  - e^x 

If  X > 0,  e^x  satisfies  (ii-a)  of  Appendix  I and  we  close  the  path  in 
the  upper  half  of  the  z-plane.  Since  e*-**  is  analytic  and  has  no  poles  in 
the  entire  finite  plane,  eq.  (I-S)  yields 


eiXx 

pClbr d* 


ilT«^  i 


X > o 


(II-l) 


Similarly  if  X < 0,  f(x,X)  satisfies  (ii-b)  and  eq.  (1-8)  is  the  one 
applicable  here.  Then  we  get 

P C 7^"  dx  " " W*5  * X < 0 ( II-2) 

When  X * 0,  f(x,X)  ■ 1,  and  although  analytic,  this  function  does  not 
satisfy  the  criteria  (ii-a)  or  (ii-b).  Therefore,  instead  of  contour 
integration,  we  go  directly  to  the  definition  (1-1)  to  compute 

pC;^"«£^o  [<1«|-*I  - 1«|-R|>  * ClnR  - lie)]  - 0 (It-3) 

e>0 


Now  we  can  combine  ( II— 1 ) , (II-2)  and  (II-3)  to  write 


Eq • (II-4)  is  valid  for  all  real  X and  £ . Separating  the  real  and  imaginary 
parts  of  (II-4),  we  obtain  two  other  results: 

P -C  ~x-r  * dx  * - tt  sgn(X  )sin(X£) 

■ "IT  sin(|X|C)  (II-6) 

P -C  *7^  dx  ■ u sgn(X  )cos(X£  ) ( II— 7 ) 

The  last  two  results  are  again  valid  for  all  real  X and  5 • If  in  (II-6) 
we  let  5 + 0 we  obtain 


p r®  cpqh*  dx  . o 
J-ro  x 


( II-8) 


This  result  can  be  also  obtained  directly  from  the  definition  of  a PV. 
Similarly,  if  we  let  £ 0 in  ( II— 7 ) we  obtain 


sinXx 

x 


dx  * wsgn(X ) 


( II-9) 


Here  the  PV  designation  has  been  deleted  since  the  integral  in  (11-9)  is 
known  to  converge  (in  fact  the  integrand  is  not  singular  at  x'a  0)  and  has 
been  evaluated  many  times  by  purely  real  methods  (Ref.  8).  This  is  not  the 
case  for  the  integrand  in  (II-8),  although  the  integrals 


r. 


S21i£  dJI,  f-b  £22ii 

X J -05  X 


are  known  to  converge  for  a,  b > 0. 

It  should  be  noted  that  making  use  of  (II-8)  and  (II-9),  we  could  obtain 
through  a change  in  the  variable  of  integration,  the  results  in  (II-6), 
(II-7)  and  ultimately  (II-4).  While  this  confirms  the  discussion  preceding 
eq.  (1-16),  it  should  be  obvious  that  the  present  procedure  is  much  faster, 
even  for  the  simple  case  considered  above. 

Plemelj's  formulas  (13)  and  (14)  yield  the  following  results: 


iX* 

f S-r~  dz 
j z-r 


-2trielX5 


X > 0 
X < 0 


-2irieiX^  H(-x) 


(11-10) 
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where 


u < 0 
u > 0 


(11-11) 


Similarly 


l\z  r 

f ^-=-  dz  - 2nielA{*  H(X) 


(11-12) 


Separation  of  real  and  imaginary  parts  in  (11-10)  and  (11-12)  yields 

/ dz  - 2irsin(X5)H(-X) 


(11-13) 


cosXz 


dz  « -2Wsin(X5)H(X) 


(11-14) 


sinXz 

z-5 


-2ncos(X5)H(-X) 


(11-15) 


sinXz 


2ttcos(X£)H(X) 


(11-16) 


A host  of  related  results  may  be  obtained  from  (II-4),  (11-10)  and 
(11-12),  if  we  admit  now  generalized  functions  as  integrands.  The  integrals 
under  consideration  are  highly  divergent,  even  when  the  singularity  at  x ■ £ 
is  absent  and  the  conventional  PV  does  not  exist.  But  we  can  extend  the  PV 
definition  to  cover  generalized  functions  as  well,  such  that  the  PV's  of 
these  integrals  will  themselves  be  generalized  functions  of  the  appropriate 
parameters . 

To  illustrate  these  ideas  let  us  start  with  eq.  (II-6)  and  differentiate 
it  once  with  respect  to  X.  By  the  preceding  discussion,  the  reversal  of  the 
P-operation  and  differentiation  in  the  l.h.s.  is  allowed  and  we  find: 


p £ H !=T  dx  * n 5x  [s«n(X)elX^] 


it  [26(X)  ♦ iCsgn(X)]  elX? 


(11-17) 


where  d/dX(sgn(X))  ■ 26 (X ) . 
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The  validity  of  interchanging  differentiation  and  the  P-operation  can  be 
verified  by  evaluating  directly  the  l.h.s.  of  ( 11-17) 


xe 


iXx 


dx  * P J_  e^x  dx  + £P 


iXx 


dx 


and  this  is  seen  to  be  identical  with  the  r.h.s.  of  (11-17). 
Proceeding  in  this  fashion  we  can  write 


k iXx 


x e 


x-£ 


dx 


(fsr)k  I”1 


isgn(X)e 


iA5 


(11-18) 


Moreover,  if  f{|(x)  is  an  N-th  order  polynomial  we  obtain  the  result 

iXx 


fN(x)e 

x-T 


dx  - f 


N (l  5")  |mi 


isgn(X  )e 


iXC 


(11-19) 


If  the  nonpolynomial  ordinary  function  f(x)  is  such  that  it  exists  as  a 
generalized  function  (Ref.  6),  then  (11-19)  can  be  extended  as  follows 


f(x)e 


iXx 


x-C 


dx  - f 


(r  a:)  1 111 


isgn(X )e 


iXf-  I 


An  example  of  a function  that  does  not  satisfy  the  above  is  f(x)  * e~^x, 
ji  > 0,  j x|  <_  « since  for  x < 0,  this  is  not  a generalized  function. 

Similar  expressions  can  be  set  down  for  (11-10)  and  (11-12)  but  we  shall 
not  do  so  here. 


2.  f(x,X)  * e^x  and  two  poles  5i  * a;  ^2  “ *> 


X i x , r i^x  „ iXx  , 


iXa  iXb 

irisgn(X)  ; a f b 


(11-20) 


The  real  and  imaginary  parts  of  (11-20)  yield 

„ i* o cosXx  N sinXa-s 

P A-  (x-a)(x-b)  dx  “ JTb 


(11-21  ) 
o 


( 11-21 x ) 
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We  assume  in  Che  above  that  aC  least  one  of  Che  two  quantities  a,b  is 


nonzero. 


The  particular  case  b ■ -a  / 0 arises  in  many  applications  and  therefore 
we  set  it  down 


p C V~2 dx  * iir£ 

x -a 


(11-22) 


Hence 


P £ dx  - -x.jnU)  iiSii 


8in( |X|a) 


(11-23) 


The  case  X ■ 0 yields,  from  eqs.  (11-22),  (11-23),  the  zero  value  for 
the  PV.  But  if  we  admit  the  possibility  that  the  arbitrary  parameter  a may 
take  the  value  a ■ 0,  we  may  ask  what  should  be  the  interpretation  of  the 
2 2 

integral  dx/(x  -a  ).  We  address  this  question  in  Appendix  VII  since  it 
affords  a good  example  of  the  6-function  interpretation  and  also  because 
this  integral  has  been  discussed  previously  in  the  literature. 

We  have  also 


-A, 


dx  ■ -2iriH(-X) 


iXa  iXb 
e - e 


(11-24) 


f ^ — 2 dx " 

-TV*  x -a 


(11-25) 


dx  * 2iriH(x) 


iXa  iXb 
e — e 


(11-26) 


; Vt  dx  - -2*H(X) 

-V*  * “« 


(11-27) 


Since  we  have  two  poles,  the  Plemelj  paths  may  be  chosen  in  four  different 
ways,  two  as  in  (11-24)  and  ( 11-26)  and  the  other  two  as  follows 
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I 


d,  - -2,1  ■i>‘l(4)  VttV< 


J 1 


e*x  . _ ...  eUaH(X)  ♦ e*bH(-X) 

■a)C^=b)  dX  2"1  Tb 


(11-28) 


(11-29) 


r «**  „ - eaaH(-A)  ♦ 

/ TT  dX  a 

w-i-»i»  x -a 


+ e'UaH(X) 


r eUx  . . .*-H<X)  ♦ e-aaH(-X) 

_W“T"2  dX  711  a 

ia/-y-*v»  x -a 


(11-30) 


(11-31) 


It  ia  of  interest  to  note  that  eqs.  (11-25)  and  (11-27)  have  no  imaginary 
parts  and  therefore  we  can  write 


; £2!ii  Jx  . 2,H(-X ) sisiA 
* i i.  a 

x -a 


; £|!i £ dx  . -2nH(X) 

* i.  l a 

-\_r>  x -a 


(11-32) 


(11-33) 


3.  f(x,A)  * e^x  and  M Poles 

Here  we  ^ant  to  evaluate  the  integral 


(11-34) 


where  g)((x)  is  an  M-th  order  polynomial  with  real  coefficients,  such  that 
all  of  its  M roots  are  simple  and  real.  If  these  roots  are  denoted  by  xm 
a well-known  formula  of  partial  fraction  decomposition  (PFD)  yields: 


i ±r-A 


1 1 


VX  m*l  g^x^)  * X® 


Then  we  inmediately  find 
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1 


(11-35) 


#iXx 


dx 


■ it  isgn(X  ) 


M 

Z 

m“l 


The  analogue  of  eq.  (11-19)  is  here 

iX  x 

dx  - f, 


— fN(x)e 

p jZ  — 


n (r  y) 


M iX  x 

* iagn(X ) Z e a 


aml  ^(xJ 


(11-36) 


If  the  roots  of  the  polynomial  g^(x)  are  not  all  real,  the  PV  can 
still  be  found,  using  the  PFD,  provided  the  locations  of  the  complex  roots 
o'f  gjj(x)  are  known.  In  this  case  contour  integration  yields  rapidly  the 
desired  result. 


Since  there  are  M poles  in  ( 11-34),  the  number  of  possible  Plemelj  paths 
is  2M,  and  hence  there  are  2M  possible  Plemelj  values  for  this  integral. 

It  is  clear  that  for  each  of  these  Plemelj  paths  ve  can  write  down  formulas 
analogous  to  (11-36). 

4.  f (x  A ) * e^x  and  an  infinity  of  poles 
The  integral  to  be  considered  here  is 


I 


r 

•LflD 


iXx 

krd* 


(11-37) 


where  g(x)  is  taken  to  be  a real  function  with  an  infinite  number  of  roots, 
all  of  which  are  real  and  simple. 

First  we  remark  that  for  X » 0,  the  PV  of  (11-37)  will  be  zero  if  g(x) 
is  an  odd  function.  In  the  opposite  case  of  even  g(x),  the  PV  may  or  may 
not  exist,  depending  on  the  detailed  behavior  of  g(x).  It  remains  to 
discuss  the  cases  X > 0 and  X < 0. 

If  we  denote  the  roots  of  g(x)  by  Xgj  we  can  evaluate  the  PV  by  contour 
integration  as  in  the  case  of  a single  pole.  Since  by  assumption  g(x)  has 
no  complex  roots  we  obtain 


..  _ iXx  ..  — _ ity  iX(x  +€e^) 

lim  p fR  *.  j_  * lim.  r . fO  £e  e m 

p 1b  I7TT  dx  + e+0  1 1 


**•  1r  gTxT 


g(« 

m 


d* 


lim  . jit  Re  e 


i0  _iXRe^ 


+ 1 J0 


g(ReI57 


d0  - 0 


(11-38) 
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This  is  Che  result  of  closing  Che  path  of  integration  in  the  upper  half 
plane,  for  X > 0. 


In  order  that  the  third  integral  vanish  in  the  limit  R*  ®,  it  will  be 
sufficient  to  have  1/|  g(z)|  ♦ a constant,  as  | z|  **’  °°,  for  all  arg  z,  such 
that  0 <_  arg  z <_  it  . In  particular  the  above  constant  may  be  zero.  This  is 
due  to  the  fact  that  in  the  numerator  we  have  exp(-XRsinP)  ♦ 0 for  R "*■  ®. 
Hence  we  obtain  from  (11-38),  using  L' Hospital's  theorem, 


iXx  iXx 


(11-39) 


If  we  assume  that  | g(z)t“^  ♦ a constant  as  | z|  + 00 , uniformly  in  arg  z, 
for  ir  < arg  z <_  2lT , we  can  extend  the  result  to  the  lower  half  plane,  to 
obtain 


eiXx 

g(x) 


dx  * -Wi 


Z 


iXx 

e m 


m g (x  ) 
m 


X < o 


(11-40) 


For  an  odd  function  g(x),  we  can  combine  these  results  into 


p£  iTU  dx 


xisgn(X) 


(11-41) 


If  the  infinite  sum  diverges,  we  may  still  interpret  it  as  a generalized 
function  of  X and  as  such  it  will  be  useful  when  appearing  as  the  kernel  of 
an  integral  operator. 

This  result  suggests  that  we  may  extend  formally  the  formula  for  the  PFD 
of  a polynomial  to  an  arbitrary  function,  even  or  odd.  To  see  this,  we 
write  for  arbitrary  g(x) 


1 

g7 x) 


-e 

m g (x  ) m 
m 


(11-42) 


where  the  roots  xm  may  be  arbitrary  complex  numbers. 

The  validity  of  (11-42)  and  its  extensions  are  discussed  in  Ref.  10. 
When  we  substitute  (11-42)  into  (11-37)  and  use  (II-4)  we  obtain  precisely 
(11-41).  Notice  that  nothing  was  said  about  the  convergence  of  the  sum  in 
(11-42). 

We  can  immediately  apply  (11-41)  to  several  functions  whose  zeros  are 
known  explicitly. 
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The  function  g(x)  ■ sin  x has  as  zeros  x„|  ■ m^;  a ■ 0,  *1*  i 
e'(x_)  - (-1)®.  As  a matter  of  fact  the  PFD  of  (sin  x)"1  is 


(11-43) 


and  g'^)  - (-1)®.  As  a matter  of  fact  the  PFD  of  (sin  x)"1  is 
well-known  (5  93  of  Ref.  9),  and  is  written  in  convergent  form  as 


i+  2 (-1)“ 

1 a-1  «2-n 


2x 

2-«V 


which  follows  directly  also  from  ( 11-42). 

The  sum  in  (11-43)  is  a generalized  function  which  can  be  written  in  a 
more  familiar  form,  when  we  use  certain  results  of  Refs.  (4)  and  (6), 


(11-44) 


A further  verification  of  (11-44)  may  be  obtained  from  the  Fourier  series 
expansion  of  the  generalized  function  (sin  x)"1.  It  can  be  shown  that 
this  (generalized)  expansion  is  given  by 


2 2 *sin(2k+l)x 

k«0 


(11-45) 


When  we  substitute  (11-45)  into  (11-37)  and  perform  the  integration  term 
by  term,  we  obtain 


(11-46) 


But  the  r.h.s.  of  (11-46)  is  identical  with  the  r.h.s.  of  (11-44).  This 
result  seems  to  imply  that  the  Fourier  transform  of  the  generalized  function 
(sinx)"1  is  identical  with  the  PV  of  the  Fourier  transform  of  the  ordinary 
function  (sinx)"1.  Since  the  following  examples  exhibit  the  same 
property,  it  may  be  conjectured  that  this  is  true  of  all  ordinary  functions 
whose  definition  can  be  extended  nontrivially  to  become  generalized 
functions. 

The  function  g(x)  ■ cosx  has  the  zeros  x^,  * (m+l/2)i*,  m ■ 0,  +1, 

+2,...,  and  g* (xjJ  ■ (-l)m+1. 


There  fore 


i\x  oo  f Is 

*Ch?, ; dx  ■ -tl-gnlx)  Z (-1)" 

m»-q» 


2u«gn(x)  2 (-1)  fi(x-(2k+l)) 

k~-a> 


(11-47) 


Now,  the  Fourier  series  expansion  of  the  generalized  function  (coax)-* 


is 


—7-  -22  (-l)k  cos(2k+l)x 
cosx  , . 

k-0 


(11-48) 


Then 


ei\x 

cosx 


dx  - 2*  2 (-1)*  [S(X+2k+l)  + <S(X-2k-l)] 

k“0 


(11-49) 


It  is  easily  seen  that  the  r.h.s.  of  (11-47)  is  identical  with  the  r.h.s.  of 
(11-49),  in  accord  with  the  conjecture  stated  above. 

The  zeros  of  the  function  g(x)  • tanx  are  Xn,  * mir,  m * 0,  +1,  +2,...  > 
and  g»  (x^,)  “ 1.  Then  we  find 


iX* 


P * cotxdx  ■ ixsgn(x)  2 e 


iXn»n 


2irisgn(X)  2 5(X-2k) 

k“-a> 


(11-50) 


Since  the  Fourier  series  expansion  of  the  generalized  function  cotx  is 


cotx  *22  sin2kx 
k"l 


(11-51) 


we  obtain 
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iXx  • 

cotx  dx  - 2iri  Z [«S(X-2k)  - «(X+2k)  ] 
k-1 


(11-52) 


which  again  is  seen  to  be  identical  with  (11-50). 

The  zeros  (of  the  function  g(x)  ■ cotx  are  Xg  ■ (m  ♦ l/2)ir,  m » 0,  +1, 
+2,...,  and  g (xm)  ■ -1.  Then 


iXx 


iX(m+4)n 


P * tanx  dx  * -irisgn(X)  Z eL*  m+2 


- -2xi  sgn(X)  Z (-1)*  6(X-2k) 

The  Fourier  series  expansion  of  the  generalized  function  tanx  is 


tanx  ■ -2  Z (-l)*1  sin2kx 
k-1 


(11-53) 


(11-54) 


and  we  find 


iXx 


tanx  dx  - -2iri  Z [S(X-2k)  - 5(X+2k)] 
k-1 

Once  more,  (11-55)  is  identical  with  ( 11-53). 


(11-55) 


All  of  the  results  (11-43)  - (11-55)  may  be  generalized,  by  analogy  with 
(11-19)  and  the  equation  following  it,  when  e^x  is  multiplied  by  a 
generalized  function  f(x),  in  the  respective  integrals. 

It  is  now  easy  to  write  the  corresponding  expressions  for  the  Plemelj 
integrals  if  we  make  use  of  the  above  results.  It  is  clear  that  there  is  an 
infinity  of  values  that  can  be  associated  with  each  integral,  since  each  of 
these  values  corresponds  to  a particular  Plemelj  path. 

We  give  one  example  where  the  function  g(x)  has  an  infinity  of  zeros, 
but  only  one  of  these  is  real,  namely  g(x)  - ex-l.  Here  the  zeros  are 
x0  - 0 (the  real  root)  and  x,,,  - 2nim,  m - +1,  + 2 . ...  . Also  we  have 
gf  (xQ)  * 1;  g^Xfl,)  - 1.  Then  the  PFD  of  (g(x)7"*  is 


x i 

e -1 


-!*;♦ 2 *,  - 


(11-56) 


m-1  x +(21011) 
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Therefore  we  obtain 


iXx 

pC  V- dx 


ex-l 


iXx  00 

- dx  + 2 I 

x m-1 


iAx 

xe 

~ 2 2 dx 

xi+(7Pmr 

(11-57) 


Ue  can  show  by  contour  integration  and  differentiation  with  respect  to  ^ 
that,  for  p > 0 


PC  T”!  dx  * iir8gn(X)e''X'p;  p>0 


(II-58) 


x +p 


where  P refers  to  the  infinite  limits  of  integration.  Using  the  above  we 
find  ■» 


P ^ — dx  * -ir6(X)  + iTTsgn(X)  +•  2nisgn(X)  I e 27f|X|m 


ex-l 


m*l 


• -tt6(X) 

- -ir6  (X ) 


+ iTTs8n(X)  [l  ♦ 2 "2] 

+ ilTsgn(X)coth(ii|X|  ) 


(11-59) 


Eq.  (11-59)  implies  also  that  for  the  integral  of  the  generalized  function 
appearing  below  we  have 


P --°^x  dx  ■ Hsgn(X)cothCn|X|  ) 
ex-l 

which  may  be  written  also  as 

sinXx  coth  y dx  ■ ifsgn(X )coth(if  | X | ) 


(11-60) 


(11-61) 


It  is  of  interest  to  note  that  the  integral  of  the  same  function  on  the 
half  range  is 

f°  -j"*  dx  • j |"coth(nX)  - ^ 1 (11-62) 

o e-l  L J 


obtained  here  through  use  of  the  PFD  (11-56)  and  of  the  result  known  from 
generalized  function  theory  that 


f*°  sinXx  dx 
J o 


(11-63) 


» 


r ; H 0 
0 ; X - 0 

The  result  (11-62)  can  be  obtained  also  by  conventional  means  since  the 
integrand  is  an  ordinary  function  without  any  singularities  (see  Ref.  9,  p. 
468).  This  again  confirms  the  legitimacy  of  the  PFD  and  its  use  under  the 
integral  sign. 

Furthermore,  we  find  from  (11-59) 

P r dx  - (X ) (11-64) 


It  is  not  difficult  to  generalise  the  results  of  (11-56)  - (11-64),  when 
g(x)  is  replaced  by  ex  - a,  0 < a,  but  this  is  left  to  the  reader. 

5.  f (x,X ) ■ e"**2;  X > 0 


The  PV  of  the  corresponding  integral  cannot  be  evaluated  by  contour 
integration.  Hence  we  use  (1-16)  to  obtain 


£ 0*  - £ e-Xx‘  sinh(2\^x)  dx 


(11-65) 


The  r.h.s.  of  (11-65)  is  obviously  an  ordinary  integral,  the  even  integrand 
has  no  singularities,  and  the  convergence  is  very  rapid. 

6.  f(x)  - e^U*2  + Xx) 


Again,  the  PV  of  the  following  integral  cannot  be  evaluated  by  contour 
integration  and  we  have  to  resort  to  other  means,  such  as  (1-16),  to 
evaluate  it.  But  first  we  rewrite  it  as 


^(yx  +Xx) 

— R 


dx 


• > 

e_l  t—  P 

4u 


,i-Ux 

x-C 


dx 


(11-66) 


where  5 * 5 + X/2y,  and  we  have  assumed  that  X and  \x  f 0,  are 

arbitrary. 


Use  of  (1-16)  in  the 
. 2 


eiyx 

x-C 


r.h.s.  of  (11-66)  yields 

dx  - ie^2  r eiMx2  sin(2y;x) 


dx 


(11-67) 


I 
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The  r.h.s.  of  (11-67)  is  a conventional  integral  with  even  integrand  that 
has  no  singularities.  Since  a closed  form  analytic  expression  for  it  is  not 
available i in  any  application  concrete  results  will  have  to  be  obtained 
either  numerically  or  by  analytic  approximation. 


The  real  and  imaginary  parts  of  C 11—67 > yield 


APPENDIX  III 


INTEGRALS  WITH  MULTIPLE  POLES 


A list  of  such  integrals  can  be  easily  constructed  from  the 
results  by  differentiation  with  respect  to  x the  function  f(x,* 
the  PV  of  the  resulting  integral,  as  explained  in  Section  III. 
shall  give  only  one  example: 


i*x 

e 

(x-S)k 


dx 


(k-1) ! 


(iM 


k-1 


■ *ik  sgn(^)^k  1 


Here  k > 0 is  an  integer. 


previous 
, and  taking 
Here  we 


(III-1) 
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APPENDIX  IV 
HALF-RANGE  INTEGRAL 

In  this  Appendix  we  present  several  examples  of  PV's  of  half -range 
integrals.  In  all  that  follows  we  assume  5 > 0. 

1.  f(x,X)  - e**»;  X / Q 

Use  of  eq.  (73)  and  (II-4)  for  F(£,X),  yields  the  expression 

_iAx  -iXx 

P dx  - iff  sgn(A  dx 

* iffsgn(X)e^  + e^  S__- — dx  (IV-1) 

Taking  real  and  imaginary  parts  we  obtain 

P dx  . -ffsin(|X|5)  - ci(|X|C)cos(XC) 

- si(|X|C)sin(|XK)  (iv-2) 

P dx  * sgn(^)  [»cosXC  + cos(XC  )si(|X|£) 

- sin(|X|Oci(|X|C)]  (IV-3) 

where  the  functions  ci(x)  and  si(x)  are  defined  in  § 8.23  of  Ref.  13. 
Incidentally  (IV-2)  and  (IV-3)  are  identical  with  the  results  given  as  No.  6 
and  No.  7 in  5 3.722  of  the  same  reference  (the  P symbol  does  not  appear 
there,  but  it  is  clear  that  the  results  are  PV  integrals). 

We  remark  that  method  (b)  of  contour  integration  is  not  applicable 
here.  The  reason  is  that  the  integral  along  the  arc  (0,  a)  increases 
without  bounds  when  | z|  ♦ •. 

When  eq.  (80)  of  method  (c)  is  used,ve  find,  after  a few  easy 
calculations 


fs-r-dx 

o x-s 


iff  sgn(X  )e 


iX5 


r 


u+iX 


du 


(IV-4) 


Since  the  PV  of  any  integral,  if  it  exists,  is  unique,  the  two  integrals  on 
the  r.h.s.  of  (IV-1)  and  (IV-4)  must  be  identical.  This  is  easily  proven, 
by  using  the  Laplace  tranform  identity  (76)  for  (x-£)~l. 


1 i 


i. 
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From  (IV-4)  we  can  also  obtain 


p r00  cosXx  m _1TSgn(\)#inX5  + J-"  — _ du 

° x_*  ° u +XZ 


(IV-5) 


• \ _ -u^ 

P /*  81  rX  dx  ■ nsgn(X)cosXC-X  J"^  ^ j du 
0 ' u +X* 


(IV-6) 


Comparing  the  results  (IV-2),  (IV-3)  with  (IV-5),  (IV-6)  respectively, 
we  notice  that  the  latter  have  monotonically  decreasing  integrands,  and  thus 
converge  much  faster  than  the  former. 

2.  f(x,X)  ■ e“^x;  X > 0 


Using  eq.  (80)  we  can  show  that  the  PV  in  this  case  is  infinite. 
Indeed,  we  obtain 

prs±4x..-*r*» 


(IV-7) 


It  is  immediately  apparent  that  the  first  integral  on  the  r.h.s.  of 
(IV-7)  is  divergent.  The  same  conclusion  can  be  obtained  from  eq.  (74), 
since  the  integral  along  the  ray  z ■ e^0^,  y > 0 and  0 < a < ff/2,  diverges. 

3.  f(x,X)  - e-**2;  X > 0 

Use  of  method  (a)  yields  here  the  result 

p J - . -2.-*2  r *-x*2  iisMiMil  dx 

J O x-c  J O X 


...  -Xx 

♦ J*°  - — f—  dx 
J o x+C 


( IV-8 ) 


We  have  used  here  (11-65). 

An  equivalent  result  may  be  also  obtained  by  method  (c),  but  (IV-8) 
seems  to  be  the  more  rapidly  convergent  expression  and  therefore  we  omit  it. 


APPENDIX  V 


# 

_ .1  _ . — 


FINITE  INTEGRALS 


In  this  Appendix  we  present  several  illustrations  of  the  results 
obtained  in  Section  V.  Here  we  assume  a < £ < b. 


1.  f(xA  ) - e*x 


Method  (a),  eq.  ( 87 )f yields 


P J 


iAx 


iAx  iA5 


>>«  dx-;be  -* 


X-£ 


dx  + e^*  log 


*5 


b-5 


2i 


f iAt  sinAt  . , b-£ 

1 J e — - — dt  + log  %— 


a± 

2 


(V-l ) 


Contour  integration  yields,  from  eq.  (88), 
iAx 


P f 


b e 


•»b+a  *\b-a  i© 

. _ iA5  . . fH  e 2 e1  2 6 ^ 

dx  - iTT  e ♦ i J Q b+a.£  -iB'  - * 

1 + — b “ e 


(V-2) 


Method  (c),  through  eq.  89,  implies  that 
iAx 


P / 


b e 


,iA5 


dx  » iir  sgn(A  )e 

r (u.iX).iXl,.~(l,-t),‘ . (u-iX).i>v(;-,)u 

“ J - 9.9 


2 ,2 

u ♦ A 


du 


(V-3) 


The  integral  on  the  r.h.s.  of  (V-3)  is  a rapidly  convergent  integral. 
Furthermore,  we  notice  that  the  first  term  on  the  r.h.s.  is  precisely  the  PV 
of  the  full-range  integral. 


Finally,  we  use  directly  the  PV  of  the  full-range  integral,  to  write 
from  eq.  (90) 


P J 


, iAx  .,r  _ iAx  _ iAx 

e dx  ■ insgn(A)el  - ££  * r~  dx  - jT  ~ir 

x-s  b X s 


x-t 


(V-4) 
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If  in  (V-4)  we  use  identity  (76),  we  can  show  that  (V-3)  end  (V-4)  ere 
indeed  identical.  But  since  the  PV  is  unique  (V-l)  and  (V-2)  oust  have  the 
same  value.  This  implies  that  we  can  establish  certain  identities 
connecting  the  various  integrals  in  the  above  four  formulas.  We  notice  that 
these  identities  are  far  from  being  trivial.  By  taking  real  and  imaginary 
parts  a number  of  other  identities  can  be  similarly  obtained. 


2.  f(x,X) 


eAX;  X real 


For  the  present  example,  methods  (c)  and  (d)  are  not  applicable, 
regardless  of  the  values  of  X.  Therefore,  we  find  from  (a)  and  (b)  the 
respective  results: 


, Xx  . XX  X5  . r t p 

p J Ihr  ■*  - / * ,'?  ■**  ♦ • l°«  £ 


2 

r Xt  sinhXt 

J e — r 

' 

2 


(V-5) 


xKi ♦ x^ie 


« Ax  < r a.  A^^“  * A « C 

« fb  e . X£  . . exp  2 2 

p ; ^ a*  • W«  % ♦ 1 J„  p b*.-2f-l6 

* 1 * 


(V-6) 


We  note  that  the  uniqueness  of  the  PV  again  will  yield  interesting 
integral  identities.  One  such  identity  is  obtained  immediately  by  equating 
to  zero  the  imaginary  part  of  the  r.h.s.  of  (V-6),  since  the  l.h.s.  is  real. 

Other  examples  may  be  treated  in  a similar  fashion,  provided  the  methods 
used  are  applicable.  That  this  is  not  always  the  case  we  have  shown  in 
example  2 above. 


APPENDIX  VI 


EVALUATION  OF  /"  dx/(x-£) 


The  integral 


(vi- O 


where  5 is  arbitrary,  is  a divergent  integral  by  the  conventional  definition 
of  Riemann  integration  as  well  as  by  its  PV  interpretation.  Moreover,  we 
note  that  the  divergence  is  not  due  to  the  presence  of  the  pole  in  the 
integrand,  since  the  integral 


VS) 


r°  dx 
a x-E 


(VI-2) 


where  |£|  < a,  is  also  divergent. 


Therefore,  the  question  arises  whether  one  can  attach  a finite 
particular  value  to  (VI-1).  Now,  we  have  shown  in  Section  iKc)  that 
(x-O"1  should  be  interpreted  as 


1 


+ C6(x-0 
P 


'32> 


If  we  perform  the  integration  indicated  in  (VI-1),  we  have 

I8^)  * IP(0  + CH(C)  (VI-3) 

To  associate  a finite  value  to  IP(£)*  we  borrow  Hadamard1  s Finite  Part 
idea  (Ref.  11)  to  calculate  it  in  the  following  fashion:  we  differentiate 
IP({[)  formally  with  respect  to  The  result  is 

4r  (IP(£)''  » J®  — 'VI-4) 

* ° (x-n2 

We  perform  the  integration  in  the  r.h.s.  of  (VI-4),  since  this  is  a 
convergent  integral.  Thus  we  obtain 

(IP(U)  « - £ (VI -5) 
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The  final  step  is  the  integration  of  the  differential  equation  (VI-5). 
Its  solution  is 


IP(5)  - B - log  Id  < 

where  B is  an  arbitrary  constant  of  integration. 

Substitution  of  (VI-6)  into  (VI-3)  yields  the  desired  result 

I8'C)  . a + CH(5)  - log  ICl  > 

But  we  notice  that,  due  to  the  properties  of  the  Heaviside  function 
H(C),  eq.  (VI-7)  can  be  written  also  as 


(VI-6) 


'VI-7) 


I8(£)  - + C2H(£)  - log  Id 


'VI-8) 


where  Ci  and  C2  are  arbitrary  constants. 

According  to  Hadamard,  what  we  have  done  here  was  to  discard  the 
singularity  present  at  upper  endpoint  of  the  interval  of  integration,  namely 


It  is  of  interest  to  note  that  we  can  obtain  result  (VI-6),  apart  from 
the  arbitrary  constant  B,  directly  from  the  definition  of  PV,  in  conjunction 
with  Hadamard* s recipe  of  discarding  the  (infinite)  contribution  from  the 
singularity  at  infinity. 


Thus,  we  can  write 


_ f40  dx  lim  lim  I rC-6  dx 

p e-0  |J-  — 


dx  . rR  dx 


log  Id  + log  |r-E 


(VI-9) 


Discarding  the  term  log  |r-C|  which  represents  the  infinite  contribution  (or 
replacing  it  with  a finite  constant  B),  we  obtain  the  result  in  (VI -6).  We 
notice  one  difference  in  the  two  procedures:  in  the  first  treatment  £ was 
arbitrary,  which  implied  that  for  Z < 0 there  was  no  pole  along  the  path  of 
integration.  In  the  PV  procedure  on  the  other  hand,  we  must  have  0 < C. 
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APPENDIX  VII 


EVALUATION  OF  £ d*/*  "5' 

In  this  Appendix  we  discuss  the  evaluation  of  the  integral  in  the 
title.  The  reason  for  doing  this  is  that  in  the  process  we  are  able  to 
emphasize  certain  points  connected  with  the  6— function  interpretation.  In 
addition,  it  appears  that  in  some  instances  the  results  given  in  the 
literature  were  erroneous.  Thus  in  Ref.  17,  this  integral  is  taken  to  be 
zero,  without  specifying  in  what  sense  or  for  what  values  of  5 is  this 
statement  valid. 


The  integral 


ICO  - £ 


dx 


2 -2 

x -5 


(VII-1) 


has  the  peculiarity  that  for  5 ■ 0,  its  PV  is  zero  as  can  be  immediately 
verified  from  the  PFD 


2 .2 
x -5 


k [;=*  ’ xic  ; 5 ^ 0 


(VII-2) 


If  5 can  take  the  value  zero,  then  we  must  write  instead  of  (VII-2)  the 
more  general  expression 


1 


2 r2 

x -5  1 g 


1 h 

2 U 


♦ K 6(5) 


+ A 5(x-5)  - ^ 

P.  P 

- B 6(x+5)l 


(VII-3) 


where  A,  B,  K are  arbitrary  constants.  Eq.  (VII-3)  follows  from  the  results 
of  Section  II  (c). 

When  we  perform  the  multiplications  indicated  in  (VII-3)  we  find 


2 -2 

x -5 


2 r2  I 

g x -5  1 P 


+ K56(5) 


__1__  I + KU^  6(e)6(x) 
x -5  1 P 


(a6(x-5)  - B6(x+5)) 


(VII-4) 
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Substituting  (VII-4)  into  (VII-1)  we  obtain: 


rg/*  \ _ -b n-  \ . - 1 


I*(g  ) - I Cg  ) ♦ C 


1 f 


C^(g) 


(VII-5) 


where  Ci  ■ (A-B)/2,  C2  * RCj  are  again  arbitrary  constants. 


If  in  (VII-5)  we  choose  lP(g  ) ■ 0,  we  have  finally  for  the  general 
determination  of  I (g  ) 


I*(g) 


:i  r 


+ C2«  (g ) 


(VII-6) 


The  variable  g in  l/g|  p takes  all  possible  values  except  g - 0. 
A particular  value  of  1(g)  is  given  in  Ref.  18,  namely 


dx 


2-2 
x -g 


*26(g) 


( VII-7 ) 


This  result,  which  corresponds  to  the  choice  Ci*0  and  in 

(VII-6),  is  not  quite  correct  in  the  way  it  is  written  in  (VII-7).  We  show 
below  in  what  sense  we  can  attach  this  particular  determination  to  1(g). 

Let  us  take  the  Fourier  transform  of  1(g): 


-iXg 


2-2 

x < 


£ dx  £ dg  £ 


-iXg 


2-2 
x -<, 


(VII-8) 


In  (VII-8)  we  have  assumed  that  inverting  the  order  of  integrations  is 
legitimate.  The  inner  integral,  over  g,  can  be  evaluated  in  the  sense  of  a 
PV  with  the  aid  of  eq.  (11-22),  which  yields  here 


-iXg 


2 -2 
x -g 


dg  - 


TTsin(IXlx) 


(VII-9) 


Substituting  (VII-9)  into  (VII-8)  we  obtain 


- » C 


Jt  sga(  | X | ) 


(VII- 10) 


I 
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Taking  now  the  inverse  Fourier  transform  of  (VII-10)  we  obtain 
IP(C)  - y £ dXelX5  *2sgn(|X| ) ■ "2  [^f  £ elX5  dx]«  *2  6(5)  (VII-11) 

All  these  manipulations  show  that  the  P symbol  attached  to  1(C)  in 
(VII-7)  applies  actually  to  a different  integral.  Moreover,  we  have  tacitly 
assumed  in  (VII-9)  that  X / 0.  The  opposite  case,  of  X « 0,  brings  us  back 
to  the  original  integral  I(x).  Thus  we  see  that  this  method  of  finding  a 
particular  determination  is  fraught  with  difficulties  and  uncertainties. 

The  conclusion  we  draw  from  the  preceding  analysis  is  that  for  any 
particular  problem,  we  should  choose  for  1(C)  (from  eq.  (VII-6))  the  value 
yielding  an  answer,  which  is  correct  according  to  criteria  imposed  by  the 
original  problem. 
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